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FOREWORD 


This  report  documents  an  in-house  effort  to  develop  more  flexible  computer 
methodology  for  optimization  or  parameterization  investigations  of  the  basic 
design  parameters  for  a  cylindrical  warhead.  The  effort  was  conducted  between 
June  and  August  1971. 
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ABSTRACT 


This  report  documents  the  Cylindrical  Warhead  Design  Optimization  segment 
of  the  Weapons  Optimization  Techniques  computer  program.  This  segment  enables 
the  user  to  optimize  or  parameterize  the  basic  design  parameters  of  a  theore¬ 
tical  warhead  for  a  given. target  or  set  of  targets.  The  warhead  lethality 
is  determined  as  a  function  of  the  basic  design  parameters:  warhead  weight, 
warhead  volume,  warhead  diameter,  charge-to-metal  ratio,  fragment  mass, 
ratio  of  warhead  length  to  diameter,  and  fragment  height-to-width  ratio. 

This  segment  can  also  optimize  or  parameterize  height  of  burst,  terminal 
velocity,  impact  angle,  and  fragment  spray  angle. 
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SECTION  I 


INTRODUCTION 

Cylindrical  Warhead  Design  Optimization  (CYDOP)  is  a  segment  of  the 
Weapon  Optimization  Techniques  (WOPS)  computer  program.  The  WOPS  r-rogram 
(Reference  1)  was  written  by  the  Martin  Marietta  Corporation  and  vas  pre¬ 
sented  to  the  Air  Force  Armament  Laboratory  (AFATL)  in  October  19  j7. 
Initially,  the  program  optimized  specific  warhead  and  delivery  parameters 
of  fragmentation,  blast,  and  fragmentation/blast  warheads  for  materiel  and 
personnel  targets  for  which  lethal  area  is  a  measure  of  effectiveness. 

Since  that  time,  WOPS  has  been  modified  several  times,  and  at  the  date  of 
this  report,  CYDOP  is  the  most  recent  of  these  modifications. 

The  first  WOPS  modification  (Reference  2)  was  a  part  of  s  follow-on 
contractual  effort.  The  contractor  modified  WOPS  to  increase  program 
flexibility  and  to  maintain  state-of-the-art  technology.  The  main  feature 
of  this  modification  is  the  capability  to  optimize  warheads  with  vai'ying 
charge-to-metal  ratios.  Other  improvements  include  the  inf orporation  of 
the  fragment  shape  (k)  factor,  the  capability  to  input  special  drag  co¬ 
efficient  tables,  the  capability  to  compute  munition  letiul  area  against 
troops  in  foxholes,  and  the  computation  of  minimum  effective  fragment 
velocities  for  personnel  and  materiel  targets  by  recent  state-of-the-art 
technology. 

In  September  1969,  a  second  modification  (Reference  5)  to  the  WOPS 
program  was  completed  by  the  contractor.  This  modififation  enables  the 
optimization  of  bomblet  pattern  size  and  the  evaluation  of  the  kill  proba¬ 
bilities  of  warheads  emitting  a  high  velocity  plane  of  fragments. 

Three  modifications,  in  addition  to  those  prepared  by  the  contractor, 
have  been  incorporated  into  the  WOPS  program  by  AFATL.  The  program  was  first 
modified  to  receive  warhead  and  target  data  fx'om  magnetic  tapes.  The  second 
modification  gave  the  capability  of  using  only  one  input  deck  to  evaluate 
similar  cases  in  the  same  computer  job. and  also  o  method  of  changing  those 
data  that  vary  from  case  to  case.  At  the  date  of  this  report,  the  third 
modification  is  under  development.  Once  completed,  the  program  will  bo 
capable  of  parameterizing  lethal  area  when  optimizing  bomblet  pattern  size 
so  that  its  sensitivity  to  lethal  area  can  bo  analyzed. 

The  WOPS  program  assumes  that  the  basic  d  s£yn  of  the  warhead  is  known 
prior  to  program  execution.  Such  warhead  cho-a*.: tfi.-tics  as  fragment  fly- 
off  velocity,  fragment  density,  and  fragment  presented  area  must  be  input 
to  the  program.  WOPS  can  assist  in  the  design  of  warheads  only  in  that  it 
can  determine  the  optimal  fragment  mass  and  the  optimal  extremes  of  the 
fragment  spray  angle  for  a  given  target  spectrum.  Often  it  is  necessary 


to  perform  analyses  to  determine  the  optimal  parameters  of  a  theoretical 
warhead.  When  a  warhead  is  in  the  design  stage,  the  characteristics  are 
not  known  but  must  be  calculated  externally  and  input  to  the  program.  The 
Cylindrical  Warhead  Data  (CYNDAT)  program  was  written  to  perform  these 
calculations . 

CYNDAT  uses  basic  design  parameters  to  compute  dimensions  and  explosive 
characteristics  of  theoretical  warheads.  Prior  to  the  CYDOP  modification, 
it  was  necessary  to  use  the  CYNDAT  output  as  the  input  to  WOPS.  Moreover, 
for  each  set  of  CYNDAT  outputs,  one  run  of  WOPS  had  to  be  made.  After  the 
CYDOP  modification,  the  final  answer  can  be  obtained  in  one  computer  run. 

This  report  documents  those  segments  of  CYDOP  that  differ  from  WOPS. 

The  documentation  consists  of  Section  II  User  Section,  Section  III  Analyst 
Section,  and  Section  IV  Test  Case.  The  User  Section  contains  a  brief 
description  of  the  program  and  an  explanation  of  the  data  input  formats. 

The  Analyst  Section  discusses  the  theory,  the  FORTRAN-source  statements,  and 
flow  charts.  Section  III  illustrates  the  input  and  gives  an  interpretation 
of  the  output  to  a  test  case.  Appendix  I  presents  the  utilization  report 
to  be  used  with  Section  II.  Appendixes  II  and  III  discuss  the  Gurney 
equation  correction  factor  and  tho  polar  zone  assumption,  respectively, 
used  in  Section  III. 


SECTION  IT 


USER  SECTION 

CYDOP  is  the  last  of  three  warhead  input  segments  to  the  WOPS  program. 

The  first  describes  the  input  for  a  warhead  with  constant  charge-to-metal 
ratio,  and  the  second  concerns  input  data  for  a  warhead  with  varying  charge- 
to-metal  ratio.  The  CYDOP  segment  uses  basic  design  parameters  to  compute 
lethal  area  so  that  the  optimal  set  of  warhead  parameters  can  be  determined. 

The  CYDOP  segment  assumes  that  the  warhead  is  cylindrical  in  shape  and 
is  primarily  a  fragmentation  bomb.  Due  to  the  methodology  employed,  the 
warhead  must  have  single-end  initiation.  Furthermore,  the  warhead  must  assume 
a  fuze-first  orientation  before  detonation. 

Assuming  that  these  conditions  exist,  the  first  step  in  executing 
the  WOPS  program  with  the  CYDOP  segment  is  to  determine  the  constraints 
defining  the  warhead.  Next,  based  on  these  constraints,  the  user  must  select 
three  independent  variables  from  the  following  set  of  warhead  parameters: 
weight,  volume,  diameter,  warhead  length-to-diameter  ratio,  and  charge-to- 
metal  ratio. 

The  user  can  then  optimize  or  parameterize  any  number  of  the  three  selected 
warhead  parameters,  the  two  fragmentation  parameters  (fragment  mass  and  frag¬ 
ment  height-to-width  ratio),  and  the  standard  warhead  parameters  (height-of- 
burst,  velocity  of  the  missile,  terminal  warhead  attack  anglo,  and  upper 
limit  of  the  first  polar  zone).  Finally,  the  output  can  bo  analyzed  to 
determine  the  optimal  sot  of  parameters  for  the  warhead. 

The  completo,  updated  version  of  the  utilization  report  for  the  WOPS 
program  is  presented  in  Appendix  I.  The  format  is  the  one  used  by  the 
Freeman  Mathematical  Laboratory  at  Uglin  Air  Force  Base.  The  program  user 
should  refer  to  this  utilization  report  while  compiling  the  data  dock  input. 

The  utilization  report  is  divided  into  Constants,  Warhead,  and  Target 
Sections.  The  Constants  Section  (Section  I)  consists  of  four  cards  that  are 
used  to  define  the  basic  parameters  describing  a  warhead -target  situation. 

The  Warhead  Section  (Section  II)  is  partitioned  into  three  subdivisions, 
only  one  of  which  may  be  used  in  a  given  run.  The  first  subdivision 
(Section  UA)  concerns  the  input  data  of  warheads  with  constant  charge-to- 
metal  ratios.  The  second  subdivision  (Section  I ID)  concerns  the  input 
data  of  warheads  with  varying  charge-to-metal  ratios.  The  third  (Section  IIC) 
describes  the  input  cards  used  in  designing  a  theoretical  warhead.  The 
Target  Section  (Section  III)  concerns  the  input  cards  that  describe  the 
target.  This  report  is  primarily  concerned  with  only  Sections  I  and  IIC. 
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The  execution  card  deck  required  to  execute  the  modified  WOPS  program 
on  the  CDC  6600  computer  is  shown  in  Figure  1.  The  data  deck  which  appears 
at  the  end  of  the  execution  deck  is  shown  in  Figure  2. 


»  Not  used  if  no  input  from  tape 


Figure  1.  Execution  Deck  for  Modified  WOPS  Program 
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*  Repeated  for  each  variable  optimized 
Section  III  deals  with  target  data  irput. 

Section  IIA  deals  with  warheads  which  have  a  constant  charge  to 
metal  ratio. 

Section  UB  deals  with  warheads  which  have  a  varying  charge  to 
metal  ratio. 

Section  IIC  deals  with  warhead  design.  The  following  is  an  illustration 
of  its  make-up: 


I 

/ . — . . . . 

Card  G 

Card  P 

r— — - - 

Card  B 

Section  IIC 

— 

Figure  2.  Data  Deck  for  Modified  WOPS  Program 
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SECTION  III 


ANALYST  SECTION 


1 .  THEORY 

In  Section  I,  it  was  pointed  out  that  a  choice  of  warhead  design  para¬ 
meters  had  to  be  made.  This  choice  is  necessary  because  there  are  five 
design  parameters  of  interest,  only  three  of  which  are  needed  to  determine 
a  solution  to  warhead  design  equations.  Therefore,  the  CYDOP  modification 
is  divided  into  10  segments  to  process  the  10  combinations  which  result  from 
five  variables  taken  three  at  a  time.  Table  I  illustrates  these  10 
combinations.  Table  II  is  a  listing  of  the  variable  names. 

TABLE  I.  COMBINATIONS  OF  THE  FIVE  DESIGN  PARAMETERS 


Combination 

Possible  Combinations  Number 


Weight 

Charge  to  metal 

Length  to  diameter 

1 

Weight 

Charge  to  octal 

Diameter 

2 

Weight 

Length  to  diameter 

Diaoetor 

5 

Weight 

length  to  diaoetor 

Volume 

4 

Weight 

Diameter 

Volume 

S 

Charge  to  metal 

Length  to  diaoetor 

Diameter 

6 

Charge  to  octal 

Length  to  diuoeter 

Volume 

7 

Charge  to  octal 

Diameter 

Volume 

3 

Volume 

Weight 

Charge  to  octal 

9 

Volume 

Diaoetor 

Length  to  diameter 

10 

The  last  two  combinations  listed  in  Table  I  do  not  uniquely  determine  a 
solution  to  warhead  design  equations;  therefore,  they  can  not  be  used. 
However,  the  first  eight  are  sufficient  to  describe  any  set  of  warhead 
constraints. 


& 


TABLE  II.  VARIABLE  NAMES 


Variable  Name 

Description 

Units 

D 

diameter  of  explosive  in  the  case 

in 

DV 

detonation  velocity  of  the  explosive 

ft/sec 

DWH 

diameter  of  warhead 

in 

FM  (1,1) 

fragment  mass 

grains 

FV  (1,1) 

fragment  fly-off  velocity 

ft/sec 

GAM4A  (or  t  ) 

Gurney  equation  correction  factor 

GC  (or  ) 

Gurney  constant  for  the  explosive 

ft/sec 

UP 

height  of  the  fragments 

in 

KDGN 

design  mode  flag 

NCOMB 

combination  number 

RCN  (or  C/M) 

ratio  of  charge  to  metal 

RHOEX 

density  of  the  explosive 

lb/  cu.  in. 

RHOSTL 

density  of  the  metal 

lb/  cu.in. 

RHW  (or  M/Hf) 

ratio  of  fragment  height  to  width 

RLO  (or  L/D) 

ratio  of  warhead  length  to  diameter 

TPRAG 

total  number  of  fragments 

TICK 

thickness  of  fragments 

in 

TKET 

constant  spray  angle 

dog 

THETAA 

computed  spray  angle 

deg 

NC 

weight  of  the  case 

lb 

WE 

weight  of  the  explosive 

lb 

WF 

width  of  the  fragments 

in 

m 

warhead  length 

in 

WVOL 

warhead  volume 

cu.  in. 

WOT 

warhead  weight 

lb 

1 

.** 


Consider,  for  example,  combination  number  1  (Table  I).  The  two  unknown 
variables,  diameter  and  volume,  are  functions  of  the  three  known  parameters. 
That  is 


DWH  =  f (WOT, C/M,  L/D) 
and 


WVOL  »  f (WOT,  C/M,  L/D). 

However,  in  combination  number  S,  neither  the  length-to-diameter  ratio  nor 
the  warhead  diameter  can  be  represented  in  terms  of  the  three  known  variables. 
That  is 


L/D  f  f (WOT,  C/M,  WVOL) 

DWH  /  f (WOT,  C/M,  WVOL). 

Similarly  for  combination  10, 

WOT  /  f (L/D,  DWH,  WVOL) 

C/M  /  f (L/D,  DWH,  WVOL). 

It  now  remains  to  derive  each  equation  for  these  combinations.  This  will 
be  done  in  combination  number  order  by  solving  for  the  two  unknown  variables 
in  terms  of  the  three  known  variables.  In  the  following  equations,  RHOSTL, 
RHOEX,  and  tt  are  constant. 


The  following  basic  equations  were  used: 


C/M  - 


WE 

WC 


WOT  -  WE  +  WC 


WVOL  a 

L/D  = 

DWH  = 
WVOL  a 


WE  ‘  WC 
RHOEX  +  RHOSTL 

NHL 

DWH 

D  +  2  (THCK) 

WHL 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 
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1)  Given:  WWT,  C/M,  L/D 
Find:  DWH,  WVOL 

From  equations  (1)  and  (2),  the  weight  of  the  case  is 


m r  -  WWT  (--t-i 

WC  '  C/M  “  1  1  (7) 

and  the  weight  of  the  explosive  (WE)  is  obtained  by  rearranging  equation 

CD. 


From  equation  (6) , 


WVOL  = 


WHL 


tt  WHL3  „ 

t  Tm2 


(8) 


So, 


WHL 


(L/D)2 

IT 


WVOL 


Combining  equation  (3)  and  (9)  gives 


(9) 


.  $ 4  (L/D)2  /  ,WE  +  WC  \ 

WHL  yRHOEX  RHOSTLy  (1° 

The  diameter  (DWH)  can  now  be  determined  by  rearranging  equation  (4),  and 
the  volume  (WVOL)  can  be  determined  from  equation  (6) . 


2)  Given:  WWT,  C/M,  DWH 
Find:  WVOL,  L/D 


The  weight  of  the  case  (WC)  and  the  weight  of  the  explosive  (WE) 
are  obtained  from  equations  (7)  and  (1),  respectively.  So,  the  warhead 
volume  (WVOL)  can  readily  be  determined  from  equation  (3), 

The  length  of  the  warhead  (WHL)  can  be  determined  by  rearranging 
equation  (6).  Then,  the  ratio  of  length  to  diameter  is  given  by  equation  (4). 

3)  Given:  WWT,  ’/D,  DWH 

Find:  WVOL,  C/M 

The  warhead  length  (WHL)  can  be  determined  by  rearranging  equation 
(4),  and  the  volume  is  given  by  equation  (6). 

The  charge- to-m*:  \  ratio  is  a  function  of  the  warhead  weight  and 
the  weight  of  the  case  (or  the  weight  of  the  explosive).  So,  the  weight 
of  the  case  must  be  determined  first.  From  equations  (3)  and  (2), 


5 


WWT-WC 


WOL  = 


wc 

RHOSTL  + 


RHOEX 


Upon  rearranging  terms, 


RHOSTL  (WVOL  •  RHOEX-WWT) 
RHOEX -RHOSTL 


The  charge-to-metal  ratio  is  then  given  by. 


(11) 


(12) 


C/M  = 


WWT-WC 

WC 


(13) 


4)  Given:  WWT,  L/0,  WVOL 
Find:  DWH,  C/M 


From  equations  (6)  and  (4), 
WVOL  « 


DWH  «  J  (DWH)3  (L/D) 


So, 


DWH  o 


V? 


WVOL 

T 7d 


Using  equation  (12),  the  charge-to-metal  ratio  can  be  determined, 


C/M 


WWT-WC 

WC 


5)  Given:  WWT,  DWH,  WVOL 
Find:  L/D,  C/M 

From  equation  (IS) , 


(14) 


(15) 


(16) 


L/D  .  £  SKS& 
tt  (DWH) 3 

Using  equation  (12),  the  charge-to-metal  ratio  can  be  determined  as  in 
equation  (13). 


(17) 


10 


6)  Given:  C/M,  L/D,  DWH 
Find:  WVOL,  WWT 

The  warhead  volume  (WVOL)  is  determined  using  equation  (14) .  From 
equations  (1)  and  (3), 


WVOL  -  WC  +  WCCC/M)  _  /  1  C/M  \ 

"  RHOSTL  RHOEX  "  iRHOSTL  RHOEXJ 

Theref oi e , 


WC 


WVOL 

1  .  C/M~ 

RHOSTL  RHOEX 


From  equations  (1)  and  (2) , 


wwt  =  wc  +  wc(c/M)  =  wc(i  +  c/m) 

7)  Given:  C/'»*,  L/D,  WVOL 
Find:  DWH,  WWT 


(18) 


(19) 


(20) 


The  diameter  is  jiven  by  equation  (15) ,  the  weight  of  the  case  is 
given  by  equation  (19) ,  and  the  warhead  weight  is  given  by  equation  (20), 

8)  Given:  C/M,  DWH,  WVOL 
Find:  _,/D,  WWT 


The  ratio  of  length  tu  dimeter  is  given  by  equation  (17).  Warhead 
weight  is  given  by  equations  (19)  and  C29). 

After  one  of  the  eight  combinations  is  executed,  the  following  warhead 
variables  are  determined:  ./eight,  volume,  diameter,  length,  charge-to-metul 
ratio,  length-tc-diameter  ratio,  weight  of  the  case,  and  weight  of  the 
explosive.  The  two  fragmentation  variables  (fragment  mass  and  fragment 
height-to-width  ratio)  are  also  determined  since  they  were  input.  It  is 
then  necessary  to  compute  the  remaining  warhead  and  fragmentation  parameters 
from  these  known  variables,  These  equations  will  be  derived  in  the  following 
paragraphs . 

Fragment  fly-off  velocity  is  determined  by  the  Gurney  formula  (derivet 
in  Reference  4).  The  equation,  for  cyiir iiical  warheads,  is 


pv  ■ m  JrzWm  (m 

where  /2E  is  a  characteristic  value  of  the  explosive,  called  the  Gurney 
constant . 
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There  has  been  much  discussion  concerning  the  accuracy  of  the  Gurney 
equation.  In  Reference  5,  a  set  of  data  points  for  correction  to  the 
Gurney  formula  is  presented  as  a  function  of  length-to- diameter  ratio.  In 
Appendix  II,  an  exponential  curve  was  fitted  to  th^ce  data  points  using 
the  least  squares  method.  The  equation  for  the  correction  factor,  y,  was 
found  to  be 


Y  =  1.0  -  (0.4486) exp [-1.2345 (L/D)]  (22) 

This  is  used  as  a  multiplier  to  the  Gurney  equation,  so  that  the  final  form 
of  the  Gurney  equation  is 


FV  -  'll  rTsCc/M)  (23) 

The  diameter  of  the  explosive,  D,  is  used  in  computing  the  spray  angle 
and  the  fragment  thickness.  The  equation  for  D  is  derived  in  the  following 
manner : 


WE  ■  Vol  ,  „  RHOEX 

(expl) 


$ 


(WHL)  RHOEX 


D 


4  (WE) 

tt  (WHL) RHOEX 


D  «  lyj- 


-w 


tt(WHL)RHOEX 


(24) 


The  thickness  of  the  fragments,  which  is  also  the  thickness  of  the  case, 
is  given  by 


THCK  =  (25) 

The  fragments  are  assumed  to  be  controlled,  rectangular  parallelepipeds. 
The  width  of  the  fragments  is  given  by 

PM  «  Vol(frag)  (RHOSTL)  (7000) 

«  THCK(WF) (HP) (RHOSTL)  (7000)  (26) 

=  THCK(WF)2(H/W)  (RHOSTL) (7000) 
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FM 

(7000) (THCK) (H/N) (RHOSTL) 

Where  7000  is  a  factor  converting  pounds  to  grains. 

The  total  number  of  fragments  is  determined  by  dividing  the  weight  of 
one  fragment  into  the  total  weight  of  all  th°  fragments, 

TFRAG  =  (7P||A3?IC,  (2; 

In  reality,  TFRAG  should  be  an  integer  since  there  should  be  an  integral 
number  of  fragments  in  a  warhead  of  given  length  and  diameter,  but  this 
fact  is  ignored  here  since  the  error  is  small. 

The  fragment  presented  area,  used  in  drag  calculations,  is  defined  as 
one-fourth  the  surface  area. 


FPA  =  ~  [2(WF)HF  +  2(WF)THCK  +  2(HF)THCK] 
=  i  [WF(HF)  +  WF  (THCK)  +  HP  (THCK] 


(29) 


Since  in  the  program  the  fragment  height,  width,  and  thickness  are  in  inches 
and  the  fragment  presented  area  is  in  square  feet,  equation  (29)  must  be 
divided  by  144.  So 

FPA  =  —•  (WP(HF)  +.WF(THCK)  +  HP  (THCK)]  (30) 


Only  ono  polar  zone  is  assumed  for  the  theoretical  warhead.  The  lower 
polar  zone  angle  is  chosen  as  90°,  or  the  angle  measured  between  the  warhead 
longitudinal  axis  and  a  line  perpendicular  to  that  axis  (Figure  3).  The 
upper  polar  zone  angle  is  computed  and,  in  general,  is  90°  plus  an  acute 
angle. 


The  upper  polar  zone  angle  is  computed  by  one  of  three  methods:  (1) 
the  angle  is  allowed  to  vary  by  optimizing  or  parameterizing  the  upper  polar 
zone  angle,  (2)  the  angle  is  computed  as  a  function  of  warhead  characteris¬ 
tics,  or  (3)  the  angle  is  input  and  held  constant.  The  first  and  third  of 
these  methods  are  self  explanatory.  The  second  method  utilizes  Shapiro's 
formula  (derived  in  Reference  5): 


THETAA  =  tan"1  j - FV(WHy - 

/ 2 (DV)  [WHL  +  (D/2)2]1' 2 


(31) 


After  computing  THETAA  by  one  of  these  three  methods,  the  upper  polar  zone 
angle,  THSUR,  can  then  be  defined  as 


THSUR  =  y  +  THETAA  (32) 

Shapiro's  formula  actually  computes  the  upper  spray  angle  of  a  single 
end  initiation  bomb  (a  in  Figure  3).  However,  the  program  makes  the 
assumption  that  representing  this  angle  by  an  angle  with  its  vertex  at  the 
center  of  the  warhead  (8  in  Figure  3)  does  not  cause  a  significant  degree 
of  error.  Appendix  III  contains  a  discussion  of  this  assumption. 


Figure  3.  Polar  Zone  Representation 
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2.  ILLUSTRATION  OF  FORTRAN  STATEMENTS 


The  WOPS  program  is  divided  into  one  mainline,  four  overlays,  and  15 
subroutines.  The  CYDOP  modification  directly  affects  only  two  overlays 
and  one  subroutine.  These  segments  will  be  discussed  (flow  charts  shown 
in  Figures  4  to  6). 

The  input  overlay  (OVERLAY  1)  is  the  first  segment  of  the  program 
affected  by  the  CYDOP  modification.  The  first  card  read  by  this  overlay 
is  in  the  input  card  type  D.  The  last  variable  read  on  this  card  is  KDGN, 
the  design  mode  flag. 

C  READ  CARD  TYPF  D 

80  READ  ( 5 i  2)  I PRNT 1 ♦ I PRNT2 » IPRNT3 ♦ I PRNT4 » I PRNT5 ♦ I DUMY* I NDXR . 

1  IDUM  »DGRID»ALPHAD»ZR*VM»R6MAX*WDI 5 ♦ KDGN 

2  FORMAT  (5I1»3I5»4F10«0*2F5*0»I5) 

C  KDGN  -  1  —  FLAG  FOR  DESIGN  MODE 

IFf IPRNT5.EQ.1 )  WRITE(6»504)  KDGN 
504  FORMAT (/♦ IX #*KDGN  =*»I5) 


If  KDGN=1,  the  program  is  to  be  executed  in  the  design  mode,  and  the  program 
assumes  that  input  Section  IIC  will  follow.  If  IPRNT5=1,  the  program 
checkout  print  indicator,  the  program  will  print  the  value  of  KDGN. 

Tn  the  next  set  of  statements  that  pertain  to  the  design  mode,  the 
first  FORTRAN  statement  determines  if  the  program  is  used  for  design.  If 
KDGN  t  1,  the  program  skips  the  design  input.  If  KDGN  =>  1,  the  program 
reads  throe  input  cards: 

. -TO-  £i<)  - . .  .  . 

C  -----  inputs  IN  OESIGN  NUDE  —  INPUT  TABLE  IIC  - - 

c  RtAO  CARD  TYPE  E 

KtA0(5,l»  WlITEt 
1  FORMAT  U2At>i 

C  READ  CARO  TYPE  F 

- -»OCtUWf  1  HE  T  t  NOOHft  - 

501  FORMAT (5F10  *0)15) 

IF<F(MO$T  .  OT.  0.)  SMOSTL  -  KMOST«12 . 096Eb 
C  REA0  CARO  TYPE  G 

HEAO<5,5ag)  ti>kK,PLO,RGM»kHH,FH(l  »1> 

502  FORMAT (7F10  ?0) 

The  first  card  contains  the  warhead  title.  The  second  contains  values  of 
parameters  that  depend  on  the  type  of  explosive  and  the  type  of  metal  used. 
Also  included  on  the  second  card  is  the  combination  number,  NC0MB,  and  the 
spray  angle,  THET,  if  it  is  held  constant.  The  third  card  contains  values 
of  design  parameters  thet  the  user  wishes  to  hold  constant  throughout  the 
program.  If  the  field  for  a  particular  variable  is  left  blank,  then  it  is 
assumed  that  the  value  of  that  variable  will  vary. 

IS 


input  Section 
IIA  or  IIB 


Input 


[•Tin 


>es 


Figure  5.  Continued 
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Figure  5.  Concluded 
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Figure  6,  Flow  Chart  -  Output 


The  next  set  of  statements  immediately  follow. and  cause  the  program  to 
print  out  the  constant  values  that  were  input. 


c  Output  constants  in  design  mode. 

WRITS (6 ,50 3 >  NCOH8.RNOEX*RHOSTUtGC*DV*TH£T.WVOL»WWT*DWM*RLD»RCM 
1  RHw.FMtUl) 

503  FORMAT ( ///*2X  **C0M8  RHO  EX  RHO  M  GURNEY  C  OET  VEL 

1  ETA  VOL  WT  OWM  R/L  C/M  H/W 

2  FM**/,lX.!S.l2F10.3t/) 


tm 
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Several  variables  have  to  be  initialized  before  continuing  the  program: 

C  S  ft  T  CONSTANT  PARAMETERS  IN  DESIGN  MODE 

NPOLAR  =  1 
MINI  =  0 
IWHDT  =  0 
NGPS(l)  =  1 
INDCD  =  2 

I F ( THET . N£ •  0.)  THE  TAA  =  THET*CONV 
ISKP  =  0 

C  END  DESIGN  MODE  INPUT  SECTION.  GO  TO  TARGET  INPUT  SECTION 

GO  TO  40n 
510  CONTINUE 


NPOLAR  =  1,  defines  the  number  of  polar  zones  as  one.  MINI  =  0  indicates 
that  the  warhead  does  net  contain  a  plane  of  fragments.  IWHDT  =  0  means 
input  Section  I IB  is  not  used.  NGPS  (1)  »  1  indicates  that  the  number  of 
classes  of  fragments  1 n  the  first(and  only}  polar  zone  is  one.  INDCD  =  2 
indicates  that  the  drag  table  for  cubes  is  used.  If  the  spray  angle  is  to 
bo  held  constant,  the  spray  angle  is  set  equel  to  the  constant  input  value. 
CONV  is  a  ccnve  ->ion  factor  to  convert  degrees  to  radians.  ISKP  =  0 
initializes  the  target  input  file  number.  GO  TO  400  causes  the  program  to 
skip  over  the  statements  used  to  input  standard  warhead  data  (input  Sections 
I  In  and  tib).  51A  CONTINUE  is  the  last  statement  in  this  segment;  it  is 
used  as  a  reference  point  if  the  design  input  is  skipped  over. 

The  next  segment  of  the  program  affected  by  the  CYDOP  modification  is  a 
computation  section,  SUBROUTINE  MOVE.  This  subroutine  is  primarily  con¬ 
cerned  with  assigning  the  new  values  of  parameters  to  proper  variables. 
However,  the  CYDOP  modification  uses  this  subroutine  to  transform  the  basic 
design  parameters  into  variables  used  in  the  lethal  area  computation.  In 
other  words,  variables  such  as  fragment  fly-off  velocity,  fragment  presented 
area,  fragment  density,  and  fragment  spray  angle  are  computed  from  the  five 
basic  warhead  and  the  two  fragmentation  parameters. 

The  assigning  of  new  values  of  parameters  to  proper  variublos  is 
controlled  by  a  DO- loop: 

DO  100  J=  1  *  I STP 
IE  (JPBK  .EQ«  2)  Isi+1 
IGO  =  NAME ( I ) 

IFUPRNT5.EQ.1  >  WRITEI6.500)  IGO 
500  FORMAT! 1X.6HIG0  =  ,13) 

GO  TO  (10,20, 30,40, 40, 40, 40*40,40, 40,40,40,40, 50, 60, 70,74* 
75,76*77  1,78, 79), (GO 

The  loop  is  executed  ISTP  times,  where  1ST?  is  equal  to  1  or  the  number  of 
variables  optimized  or  parameterized.  The  value  of  ISTP  is  determined 
depending  on  conditions  set  in  the  mainline.  NAME  (I)  contains  a  value 
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representing  a  variable  that  is  optimized  or  parameterized.  Depending  on 
this  value,  the  program  branches  to  one  of  13  different  segments.  State¬ 
ments  10,  20,  30,  50,  and  70  deal  only  with  standard  variables  (namely, 
HOB,  VM,  ALPHA,  ZONEL,  and  MlWI-mass).  If  KDGN=1,  statements  40  and  60 
transfer  immediately  to  statements  73  and  72,  respectively.  This  is  done 
because  these  sections  are  not  used  if  the  CYDOP  segment  is  executed. 


AO  CONTINUE 

IF  (KDGN»FQ'>1  )  GO  TO  70 

•  •  • 

•  •  » 

•  •  * 

60  CONTINUE 

IF(KDGN*EQ.1>  GO  TO  72 

•  •  • 

%  •  • 

t  9  ♦ 

72  TH5U( 1 l  =  X < I > 

UiSOR(l)  =  XU)*CONV 
GO  TO  100 

73  PM<  1,1)  =  X(  I) 

GO  TO  100 


Statement  numbers  74  through  79  assign  design  parameters  to  proper 
variables,  Statement  100  is  tho  end  of  the  DO -loop  that  controls  the 
switching  of  variables. 

7 A  WVOL  -  X( n 
GO  TO  100 

75  WWTsXI  I  ) 

GO  TO  100 

76  OWH  =  XII) 

GO  TO  100 

77  RID  6=  XU  ) 

GO  TO  100 

7fl  RCM  a  XU) 

GO  TO  100 
79  RHW  a  XII) 

100  CONTINUE 

The  next  two  statements  cause  tho  program  to  skip  tho  design  computation, 
if  it  is  not  applicable. 


C  SKIP  DESIGN  COMPUTATION  IF  NOT  APPLICABLE' 

It (KDGN.NE.l )  GO  TO  110 

1F<  JP0K.NE.2^ND.  J  .NS.  1. ANB.NPAR.GT »0 )  GO  TO  110 
8  •  RHOEX 

IFUPRNT5-E0.1 )  WRITE(6»501)  NCOMB 
501  FORMAT (1X**NC0MB 
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B  is  set  equal  to  RHOEX  for  ease  in  reading  the  equations  that  follow. 

If  IPRNT5=1,  the  value  of  the  combination  number  is  printed. 

Immediately  following  the  previous  statements,  the  program  branches 
into  one  of  ten  segments  depending  on  the  combination  number.  These 
segments  correspond  to  the  combinations  described  in  Table  I. 

0  GO  TO  SEGMENT  UfcFIWEO  GV  THE  COMBINATION  NUMBER, 

o  IN  EACH  0  LMBI N A  1 1  ON  NUMBER  SEGMENT,  DETERMINE  THE  V  ALUE 

c  for  THE  TWO  UNKNOWN  VARIABLES  IN  TERMS  0F3  KNOWN  VARIABLES 
GOTO (8i, 82, 83,*W,  85 ,8b  ,  87 , 88 , 8 9, 89>  ,NCOMB 


Combination  number  1  determines  the  diameter  of  the  warhead  and  the 
volume  of  the  warhead  in  terms  of  the  charge-to-metal  ratio,  length-to- 
diameter  ratio,  and  warhead  weight: 


C  COMBINATION  1 

81  WC  =  WWT / I RCM+1 •  ) 

WE  =  WC*RCM 

WHL  =  ( ; (WE/BI+I WC/RHOSTL))*(RLD**2 )#4*/PI )**ONETRD 

DWH  «  WHL/RLD 

WVOL  =  (PI*DWH**2 J*WHL/4. 

GO  TO  93 


The  weight  of  the  case,  the  weight  of  the  explosive,  and  the  warhead  length 
are  determined  as  intermediary  steps,  so  they  need  not  be  determined  in 
lator  stops. 

Similar  statements  are  used  for  the  remaining  combination  numbers, 
solving  for  the  two  unknown  parameters  in  terms  of  the  three  known  variables. 
If  NC0M8  is  input  as  9  or  10,  an  error  code  is  printed  and  the  run  terminated. 

0  COMBINATION  2 

62  WO  =  W«T/(KGN*i.) 

HE  a  WC’kCM 

WVUL  -  (wC/KfOSTl) 

WHL  a  4,*WVOL/(RI*eHH**2) 

RlU  -  HhL/UWH 

G ;  TO  93 

C  COMBINATION  3 

63  WHl  *  v  .O*OHf 

WVOL-  \  '4H**2>#WriL*P 1/4* 

-  We“  -  -9*t<KHOSTt^N¥i)L-*tWT  )/(RMOSTt-tiM - -  - - - 

ROM-  ( Wh  T“WO ) /WO 

.  SO  TO  92 _  _ _ _ _  _ 
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C  '  CORBIN  AT 1UN  4 

84  DWn  =  <4.*WV0L/ (P1*RLU>  )**GNETRD 

WC  =  B* ( (KHbSTL*WVOL-WMT)/ (RHOSTl-B) ) 

- - RCH*  <mnt-wc-)vnc- - -  - - 

GO  TO  91 

C  COMBINATION  5 

65  RLU  =  4.*WVCL/(PI*0WH**3) 

WC  *  9*T{RH0STL*WV0L-NWT)/(RHGSTL-8>  > 

RCM  =  (WWT-WU/WC 

-  . . 90  -TO  91 - - - -  -  —  - - - 

G  COMBINATION  6 

86  WVOL  =  PI*RLO*<OWH**3>/4. 

WC  =  WVOL/ (  (RCM/BH- (l./RHOSTL)  > 

WWT  *  WC* (ROM+i t )  . 

GO  TO  91 

€ . . COMBINATION-  ?- . -  -  - 

8 7  QWH  =  (4.*WV0t/<Pi*RLD) >**ONETRD 

WWT  a  (WV0L/<<RCM/6>  +  (i./RHQSTLn  >*{1.+RCM> 

GO  TO  90 

C  COMBINATION  8 

88  RLO  =  4.*WVUL/(Pl*0NH**3) 

GO  TO  90 

C  COMBINATIONS  9  ANC  10 

C  T Hii  UNKNOWNS  IN  COMBINATIONS  9  AND  10  CANNOT  BE  REPRESENTED 

0  IN  TERMS  OE  THfc  KNOWN  VARIABLES 

69  WRITE(6,i9> 

-  ^-^OMNMTM-lR9rTAJUt^W-0-~&feT>fcNO€NT  VA-RlABifcS  INPUT-IN  OE SIGN- MODE** *9 

CALL  EXIT 

In  some  of  these  combinations,  weight  of  the  case,  weight  of  the  explosive, 
and  warhead  length  are  not  determined  as  intermediary  steps.  The  following 
statements  are  included  to  compute  them  when  necessary: 


C  DETERMINE  W/H  CHARACTERISTIC  VALUES  CALCULATED  IN  SOME  OF 

C  THE  ABOVE  SEGMENTS  BUT  NOT  CALCULATED  IN  OTHERS 

90  WC  «  WWT / { RCM+1 • ) 

91  WHL=  RID*DWH 

92  WE  =  WWT-WC 


The  program  must  now  transform  these  design  parameters  into  variables 
used  in  the  lethal  area  computation.  The  first  variable  computed  is  the 
fragment  fly-off  velocity.  The  equation  utilizes  the  Gurney  formula  with 
the  exponential  correction  factor  derived  in  Appendix  II: 

C  DETERMINE  REMAINING  W/H  CHARACTERISTICS 

93  GAMMA  =  1*0  -  *4486*EXP( -1 *2345*RLD ) 

FV<1*1J  »  GAMMA*GC*SQRT(RCM/( l*+*5*RCM) ) 
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The  variable  D,  diameter  of  the  explosive,  is  computed  next. 

0  =  SGRT{DWH**2-{4.*WC/(PI*WHL*RH0STl> > ) 

THCK  =  ( DWH-D ) / 2* 

WF  =  SQRT< FM( 1* 1 )/< 7000**THCK*RHW*RHOSTL) ) 

HF  a  WF*RHW 

TFRAG  =  7000.*WC/FM( 1 ,1) 


The  fragment  thickness,  width,  and  height  are  then  determined.  TFRAG,  the 
total  number  of  fragments,  is  determined  by  dividing  the  weight  of  one  fragment 
into  the  total  weight  of  the  case. 

Since  a  new  warhead  length  is  computed  each  time,  the  correction  factor 
to  the  height  of  burst  also  changes. 

C  HOB  =  HOB  +  PRESENT  CORRECTION  FACTOR  -  OLD  CORRECTION  FACT 
ZB  =  ZB  +  (SALP*WHL/24. J  -  SALP*WDI$F 
WDISF=  WHL/24. 

The  height  of  burst  must  be  adjusted  accordingly.  This  is  done  by 
subtracting  the  old  correction  factor  and  adding  the  new  factor  to  the 
height  of  burst.  The  distance  from  the  center  of  the  warhead  to  the  nose 
of  the  warhead,  WDISF,  must  also  be  corrected.  This  distance  is  one-half 
the  warhead  length  (divided  by  12  to  change  inches  to  feet). 

The  fragment  presented  area  is  computed  by  raking  one-fourth  the  surface 
of  the  fragment.  This  is  then  divided  by  144  to  convert  square  inches  to 
square  feet. 

C  FPA  IS  ONE-FOURTH  SURFACE  AREA  DIVIDED  BY  144—  TO  CHANGE  SQ 

C  IN  TO  SQ  FT 

FPA (1*11  =  ( WF *HF+WF* THCK +HF# THCK ) /288 • 

The  upper  polar  zone  angle  is  determined  by  one  of  the  three  methods 
described  earlier.  If  IZN=1,  the  upper  polar  zone  angle  is  optimized  or 
parameterized.  If  THET=0,  the  upper  polar  zone  angle  is  held  constant. 

If  neither  of  the  above  two  combinations  exist,  the  spray  angle  is  computed 
by  Shapiro's  formula  and  defined  as  the  upper  polar  zone  angle. 

THREE  OPTIONS  TO  COMPUTE  SPRAY— 

1)  IF  UPPER  ZONE  IS  OPTIMIZED  UZN*1)  THEN  THr T A  IS  DETERMINED 
BY  PRESENT  VALUE  OF  THSUR I  UPPER  POLAR  ZONE  IN  RADIANS) 
IFUZN.EQ.l)  GO  TO  95 
C  2)  THETA  IS  INPUT  AND  HELD  CONSTANT  <  THET  .NE.O*) 

IF(THET.NE,0. )  GO  TO  94 

C  3)  THETA  IS  COMPUTED  BY  SHAPIRO  S  FORMULA 

THFTAA  =  ATANtFV<l*l)*WHL/<2**OV*SGRT(WHL**2+<D/2.)**2) )) 

94  tFlTHET.NE.O*)  THETAA  =  .TH£T*CONV 


The  next  set  of  statements  define  the  polar  zone  in  terms  of  the  variable 
used  in  the  lethal  area  computation.  The  last  statement  in  this  set  converts 
the  spray  angle  to  degrees. 

THSUR(l)  =  HFP I  +  THETAA 
THSU(l)  =  THSUR( 1 )/CONV 
95  THSLR(l)  «  HFP I 

THSL(l)  =  HFPI/CONV 

IFU2N.EQ.1)  THETAA  =  THSUR  (  1  ) -HFP  I 

THETAA  =  THETAA/CONV 

After  the  density  of  the  fragments  is  determined,  the  subroutine  returns 
control  to  the  mainline. 


DENOM  =  P I 2*( COS ( THSLR ( 1 ) ) -COS( THSUR (111) 
RHO(l,l)  =  TFRAG/DENOM 
110  RETURN 
END 


The  design  mode  output  is  executed  in  program  LETHAR  (0VERLAY2) .  This 
overlay  computes  the  probability  of  kill,  prints  the  effective  and  lethal 
area  table,  and  prints  the  PK  matrix.  The  design  output  is  placed  immediately 
after  the  effective  and  lethal  area  table  (if  the  effective  and  lethal  area 
output -is  requested).  If  none  of  the  print  indicators  are  set,  the  inputs, 
the  design  mode  outputs,  and  the  summary  will  be  the  only  items  printed. 


C  PRINT  EFFECTIVE  AND  LETHAL  AREA  DATA 

1605  WRITE  (6*111 

WRITE  (6*12)  (  PLEVEL(M) *XEA( M* I ) *XLA (M* I ) *XPK(M* I  1 *M  =  1*10) 
WRITE  (6*13)  PLFVEL( 11 ) »XEA( 11 *1) ,XLA( 11* I ) ,XPK( 11* I ) 

WRITE  ( 6  ♦  1 A 1  PLEVEL ( 12 1 »XEA( 12* I  1 *XLA( 12* I ) »XPK ( 12  ♦  I) 

WRITE  (6*12)  PLEVEL ( 13 ) *XEA( 1 3* I ) *XLA( 1 3* I ) *XPK( 1 3*  1 1 
C  PRINT  DESIGN  MODE  COMPUTATION 

1610  IF(KDGN.FO.I) 

TH r„  XWRITEC6.101)  NTR  * WWT ♦ WVOL * RLD * RCM* RHW ♦ FM ( 1*1 1 *WE*WC*WHL.DWH* 
tHCK  1 *HF,WF.TFRAG.FV( 1*1 ), THETAA 

101  FORMAT t /*4X»*NUM  WWT  •  WVOL  L/D  C/M  H/W  FM 

1  WE  WC  LWH  DWH  FTK  FHT  FWD  FNUMB 

2  FVEL  ANGL»/4X*I3*2F8.2»3F8.3*3F8.1,2F9.3*3F8.3*2F8«0*F7.1) 


SECTION  IV 


TEST  CASE 

Table  III  illustrates  the  input  deck  for  a  typical  execution  of  the 
CYDOP  computer  program.  Five  variables  were  parameterized,  and  a  total 
of  162  warheads  were  evaluated  against  a  single  target  element  in  84  seconds 
of  central  processor  unit  time  on  a  CDC  6600  computer. 

Table  IV  illustrates  a  portion  of  the  output  obtained  from  this  test 
case.  This  output  is  optional  and  is  obtained  by  setting  IPP,2  equal  to  1 
on  card  D.  It  is  strongly  recommended  that,  when  using  the  warhead  mode, 
this  output  be  obtained  since  it  contains  information  about  the  warhead 
that  will  not  be  found  in  the  summary.  The  summary  is  always  output  and 
is  illustrated  in  Table  V. 

There  are  several  variables  in  Table  IV  which  have  not  been  encountered 
previously.  Effective  area  is  the  amount  of  area  in  which  there  is  a  PK 
equal  to  or  greater  than  the  designated  PK  (Minimal).  PK  (Mean)  is  defined 
as  follows: 

//PK  da 

PK  (Mean)  = 

Effective  Data 

where  EA  is  the  Effective  Area'.  Lethal  Area,  LA,  is 
LA  =  PK(Mean) *EA 

In  the  final  line  of  the  output,  there  are  several  other  new  variables: 

a.  NUM  is  the  number  with  0  indicating  the  first  case. 

b.  FTK.  is  the  fragment  thickness, 

c.  FHT  is  the  fragment  height. 

d.  FWD  is  the  fragment  width. 

e.  FNUMB  is  the  number  of  fragments. 

f.  FVEL  is  the  initial  velocity  of  the  fragments. 

g.  ANGL  is  the  upper  angular  limit  of  the  fragment  spray. 

Table  V  contains  the  summary  output  for  the  first  12  cases.  Not  all  cf 
the  variables  being  parameterized  are  shown  in  the  summary.  This  deficiency 
is  inconvenient,  and  the  program  is  being  modified  to  correct  the  problem. 


28 


351255V1# 


TABLE  IV.  TEST  CASE  FOR  DESIGN  MODE 

TESI  CASE  WARHEAD  TARGET  ELEMENT  NO.  1 

WARHEAD  ENTRY  ANGLE  6D.0  DEGREES 

WARHEAD  VELOCITY  250.0  FT/SEC 

FM  1  15.0 

DURST  HEIGHT  .198  F cET  - — ' ~~ . * 


PK 

(MINIMAL) 

EFFECTIVE  AREA 
<sq  FT) 

LtTH Al  AREA 
(SQ  FT) 

PK 

(MEAN) 

1.0  - 

OVO0000--  "■ 

0  «  0000-0’- 

■  -  - . . 

tr.  0000-00-0  0 

.9 

352.33701 

338.09014 

.95956465 

.8 

352.33701 

338.09014 

.55956465 

.7 

352.3-3701 . 

330T000I4  - 

- — 

.959564-65 

.6 

352.33701 

338,09014 

»  55956h65 

♦  5 

352.33701 

338.09014 

.95956465 

'  .4 

“  352.33701 - 

330.-09014  - 

— — - - - 

■’  .  559504  65 

.3 

352.33,701 

*  338.09014 

. 9595b465 

.2 

827.70571 

465.05914 

.56186533 

.1  '  - 

1499*7366 1‘ - 

‘  "575;'45439 

- ..  -  .. . 

.  38370364' 

.005 

2429,06160 

662.59808 

.  2/27  7945 

.  0001 

2429,06160 

662.59808 

.  27  27  7545 

0.0 

55115.62613 

662.59008 

- 

•  01202196 

NUM  WWT 

WVOL  L/O 

C/M  H/H 

FM” 

u  -  T.Qtt  t.50  1,000- - 

‘  .  300  1.000 

*  “15V0 

.2' 

WC  LHH 

DWH  FTK 

FMT  FWD 

FNUHB 

FVEL  ANGL 

‘"*8  “'"2.“023‘  • 

— 2,023-  “rttrO - 

•  IT* - ~.170— 

-  35"9i 

3900.  -  4.1 
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TABLE  V.  TEST  CASE  FOR  DESIGN  MODE 


SUMMARY 


TARGET 

NO.  1 

TARGET 

p 5  V4 

WOOD 

TRIAL 

TARGET 

ALPHA 

HOB 

W/H 

Fid) 

LETHAL 

NO. 

NO. 

(DEG) 

(FT) 

V  EL. 

GRAINS 

AREA 

(FPS) 

SQ.FT. 

0 

1 

60.  C  ' 

- ,  2- 

“■250.13 

15.0 

55T. 

i 

1 

75.0 

.2 

25  0 . 0 

15.0 

l<to4. 

2 

1 

90.0 

.  2 

250.  0 

15. 0 

0. 

3 

1 

60.0 

.2 

25  0. 0 

30.0 

60  2. 

4 

1 

75.  0 

.2 

250.0 

30.0 

135  0* 

5 

1 

90.  0 

.2 

250.  0 

30.0 

0. 

6 

~i . 

6  d .  C 

’  ".2 

250.0 

15. 0 

940. 

7 

l 

75.0 

.2 

25  0.  0 

15.0 

2011. 

8 

l 

90.0 

.2 

25  0.  0 

15. 0 

0. 

9 

l 

bO.  0 

.2 

250.0 

30.  0 

862. 

10 

l 

75.  0 

.2 

250,0 

30.0 

1900. 

11 

l 

90.0 

.2 

250.  0 

30.0 

0. 

However,  by  simply  referring  to  Card  C  of  the  input  deck,  the  values  of  the 
variables  being  parameterized  can  be  determined  for  any  of  the  cases.  The 
rule  is  that  the  variable  input  on  the  first  of  the  Type  C  cards  will  change 
most  quickly  while  the  variable  input  on  the  last  of  the  Type  C  cards  will 
change  most  slowly.  In  this  test  case,  therefore,  ALPHA  varies  most  quickly 
and  L/D  changes  most  slowly. 
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SECTION  V 


CONCLUSIONS  AND  RECOMMENDATIONS 


1 .  CONCLUSIONS 

a.  The  Weapon  Optimization  Techniques  program  has  been  modified  so 
that  it  is  possible  to  optimize  or  parameterize  basic  design  parameters 
for  cylindrical  warheads. 

b.  The  methodology  for  fragment  fly-off  velocity  has  been  investigated, 
and  a  correction-factor  equation  to  the  Gurney  formula  has  been  introduced. 

2 .  RECOMMENDATIONS 

a.  The  possibility  of  considering  blast  effects  when  using  this  program 
to  optimize  design  should  be  investigated. 

b.  A  design  option  to  this  program  for  alternative  warhead  shapes 
should  be  considered  as  a  possible  future  modification. 

c.  The  possibility  of  expanding  the  warhead/ target  representation  to 
consider  more  than  just  point  locations  should  be  investigated. 

d.  Future  programming  efforts  should  be  thoroughly  documented  to 
eliminate  any  uncertainty  as  to  how  the  program  accomplishes  its  computations. 
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SECTION  I 


constants,  section 


CARO  A 


column  item  scaling 

1-72  GTITLE  «  GENERAL  TITLE  12A6 

CARO  8 

COLUMN  ITEM  SCALING 

1-5  NVOPT  ®  NUMBER  OF  VARIABLES  TO  BE  OPT IM1ZEOIMAX  OF  15)  15 

6-10  NPAR  =  NUMBER  OF  VARIABLES  TO  BE  PARAME TER i ZED ( MAX  OF  8)  15 

IF  NVQPT  IS  GREATER  THAN  0*  NPAR  »  0  AND  VICE-VERSA. 
11-20  CC1  F10.0 

21-30  CC2  F  1 0 • 0 

OPT  HI  2  AT  I  ON  PROCEDURE  WILL  NOT  CONVERGE  UNTIL  THE 


FOLLOWING  TWO  CONDITIONS  ARE  SATISFIED.  AL t  I  )  IS  LETHAL 
AREA  OF  TRIAL  I. 

1.  ABSOLUTE  VALUE  OF M At ( N- 1 > -AL I N-2 > ) /AL < N-2 ) ) 

MUST  BE  LESS  ThAN  CCi  AND 

2.  ABSOLUTE  VALUE  OF!  CALIN)-AHN-l)  l/ALlN-l)  ) 

MUST  BE  LESS  THAN  CC2. 

35  I  PR 6  *  It  PRINTS  INTERMEDIATE  VALUES  IN  OPTIMISATION  11 

PORTION  OF  PROGRAM. 

AO  tQPPAT*  l,  OPTIMIZES  PATTERN  SIZE*  II 

AS  I  EACH  a  1,  IN  PARAMETRIC  MODE. INDICATES  PATTERN  SIZE  WILL  11 

BE  OPTIMIZED  FOR  EACH  SET  OF  CONDITIONS. 

0.  PATTERN  SUE  WILL  OE  OPTIMIZED  FOR  THE  LARGEST 
VALUE  OF  N 

PI  I AL 1  WHERE  N  IS  THE  NO.  OF 

I  target  elements. 

46  IAPAR  ■  It  lETHAL  AREA  WILL  BE  PARAMETtRUED  II 

NOTE.  IF  LAPAR  *  1.  SKIP  TO  CARD  Wl* 


CARO  C 


COLUMN  ITEM  SCALING 

1-6  NAME  I  I J  •NAME  Of  VARIABLE  TO  BE  OPTIMIZED  OR  A6 

parameterized  from  variable  TABLE 

7~r<  Blank  <.x 

11-20  XMNII)  a  VARIABLE  minimum  F10.0 

1-30  XMKClI  *  VARIABLE  MAXIMUM  FjQtO 

1-40  xin  =  FIRST  ESTIMATE  IN  OPTIMIZATION  MODE  F10.0 

«  DELTA  X  IN  PARAMg TR SC  MODE 


CARO  C  IS  REPEATED  FOR  EACH  VAR! ABLE  BEING  OPTIMIZE* 
OR  PARAMETERIZED 
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VARIABLE  TABLE 


COLUMN 

ITEM 

dimension 

1-3 

HOB  = 

WARHEAD  HEIGHT  OF  BURST 

FT. 

1-2 

VM  = 

TERMINAL  WARHEAD  VELOCITY 

FT. /SEC. 

1-3 

ALPHA 

=  terminal  warhead  attack  ANGLE 

DEGREES 

1-5 

FM  (  I  ) 

=  FRAGMENT  MASS,  GROUP  I  (MAX  OF  10  GROUPS) 

GRAINS 

1-5 

ZONEL 

=  LOWER  LIMIT  OF  FIRST  ZONE  IN  WARHEAD 

DEGREES 

1-5 

ZONEU 

=  UPPER  LIMIT  OF  LAST  ZONE  IN  WARHEAD 

DEGREES 

1-3 

FM1 

=  FRAGMENT  MASS  FOR  PLANE  OF  FRAGMENTS 

GRAINS 

1-3 

* 

VOL  = 

VOLUME  OF  WARHEAD 

CU.  IN. 

1-3 

* 

WWT  " 

WEIGHT  OF  WARHEAD 

LBS 

1-3 

* 

DWH  = 

DIAMETER  OF  WARHEAD 

IN. 

1-3 

L/3  = 

RATIO  OF  LENGTH  TO  DIAMETER 

- 

1-3 

* 

?./'■*  - 

RATIO  OF  WEIGHT  OF  CHARGE  TO  WEIGHT  CF  r.ETAL 

- 

1-3 

* 

H/W  = 

RATIO  OF  HEIGHT  TO  WIDTH  OF  FRAGMENTS 

- 

NAMES  MUST  BEGIN  IN  COLUMN  1  AND  BE  IN  THE  EXACT  FORMAT 
AS  INDICATED  IN  THIS  table 

*  USE  WITH  DESIGN  MODE  ONLY  ( KDGN  =  1  ).  DESIGN  MODE  CAN  ALoO 
USE  HOB,  VM»  ALPHA,  FM ( 1 ) ,  AND  ZONEU.  REFER  TO  NOTE  AT  I  HE 
BOTTOM  OF  THE  FIRST  PAGE  OF  SECTION  I  I C  BEFORE  CHOOSING  ANY 
OF  THE  DESIGN  PARAMETERS. 


CARD  D 


COLUMN 

1 

2 

3 

4 

5 

IS 


ITEM  SCALING 

IPRNT1  =  1,  PRINTS  STATIC,  DYNAMIC  AND  FRAGMENTATION  II 

WARHEAD  DATA 

IPRNT2  =  1,  PRINTS  EFFECTIVE  AND  LETHAL  AREA  DATA  U 

IPRNT3  =  1,  PRINTS  PK  MATRIX  '  II 

I  PR, NT  4  =  1,  PRINTS  TARGET  DATA  INPUT  11 

IPRNT5  =  1,  SPECIAL  PRINT  USED  FOR  PROGRAM  CHECKOUT  II 

SET  THE  INDICATOR  TO  0  TO  OMIT  THE  SPECIFIED  PRINTOUT 

INDXR  =  1,  INDICATES  EVERY  RADIUS  IN  THE  TABLE  OF  RADII  IS 

FOR  THE  TARGET  DAMAGE  ASSESSMENT  uRJD  IS  TO  BE 


USED 

=  2,  INDICATES  EVERY  OTHER  RADIUS 
=  3,  INDICATES  EVF-RY  THIRD  RADIUS 


ETC 

21-30  DGilP  ■=  ANGULAR  INCREMENT  FOR  TARGET  DAMAGE  ASSESSMENT  F10,0 


GRID.  SHOULD  DWIDE  INTO  1B0  EVENLY  (DEGREES) 

31-40  ALPHA  =  WARHEAD  TERMINAL  ATTACK  ANGLE  (DEGREES)  FIO.O 

41-50  H03=  WARHEAD  HEIGHT  OF  BURST  (FT)  FIO.O 

51-60  VM  -  WARHEAD  TERMINAL  VELOCITY  (FT/SEC)  FIO.O 

61-65  RGMAX  =  MAXIMUM  EFFECTIVE  RADIUS  TO  BE  EVALUATED.  IF  F5.0 

RGMAX  =  0  .  GRID  RADIUS  WILL  BE  Dt  TERM! NED  BY  A 
CALCULATED  EFFECTIVE  RANGE  OF  WARHEAD. 

66-70  WDIS  »  DISTANCE  FROM  ThE  NOSE  OF  THE  MUillt  TO  THE.  FS.O 

CENTER  OF  THE  WH  SECT  I  ON 1 1 NCHES ) 

75  KDGN  =  1,  DESIGN  MODE  FLAG.  USE  SECTION  11  C  IF  KDGN = 1 •  15 
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SECTION'  1 1 


WARHEAD  SECTION 


COLUMN 

1 

2-3 

4-9 


COLUMN 

1-72 


COLUMN 

1-5 

6-10 


11-20 

U-30 

35 

36-60 

61-70 

71 


72 


NOTE-  USE  TABLE  IIA.IIB*OR  1  I C  AS  APPROPRIATE 


SECTION  I  I  A 

(WARHEADS  WITH  CONSTANT  CHARGE  TO  MASS  RATIO) 


CARD  t 


ITEM  SCALING 

ITAPt  «  1.  INDICATES  CERTAIN  VARIABLES  WILL  BE  READ  11 

FROM  TAPE 

ISKP  =  NO.  OF  THE  FILE  TO  BE  READ  FROM  TAPE.  12 

TITLE  =  ABBREVIATED  NAME  OF  THE  FILE  TO  BE  READ  FROM  A6 

TAPE. 


CARD  F  -  OMIT  WHEN  ITAPE=1 

ITEM  SCALING 

WTITLE  =  WARHEAD  TITLE  12A6 


CARD  G  -  WHEN  ITAPEbU  INPUT  COLUMNS  11-35  ONLY 


ITEM  SCALING 

BLANK  IS 

NPOLAft  »  NUMBER  OF  POLAR  ZONES  REQUIRED  TO  DESCRIBE  THE  15 

FRAGMENTATION  DATA  -OR  WARHEAD  (DEG)  (MAX  OF  36) 


NOTE.  IF  WH  CONTAINS  .ONLY  A  PLANE  OF  FRAGMENTS 

ANO  NO  FRAGMENTS  IN  POLAR  ZONES*  THEN  NPOtAR 
MAY  BE  SET  TO  0  AND  CARDS  H  AND  1  OMITTED. 

PBEXP  a  THE  BARE  CHARGE  OF  TNT  EQUIVALENT  EXPLOSIVE  IN  F10.0 
THE  WARHEAD  (L»S) 

rhomtl  °  density  of  fragment  MATERIAL  U8S/!N**3>  F 10.0 

(NOT  USED  IF  FRAGMENT  IS  STEEL) 

MINI  »  1*  INDICATES  WH  CONTAINS  A  PLANE  OF  FRAGMENTS.  15 

blank 

FACTK  ®  K  FACTOR*  ( FT »• 2 /6R**Z/3 )  F10.0 

(FRAG  FRAGMENT  OENSlfY  INDICATOR  II 

*  0  the  fragment  density,  RHO ( 1 1 J ) »  *S  IkPUT  in 
FRAGMENTS/STEP AOt An 


»  1  THE  FRAGMENT  DENSITY  .'5  INPUT  AS  THE  TOTAL  NUMBER 

of  fragments  and  the  program  computes  fragments/ 
steradian 

UDCO  >*  FRAGMENT  DRAG  indicator,  1  ■  SPHERE.  2  »  CUBE.  II 

3  a  RANDOM,  IF  UDCO  GT  3,  A  SPECIAL  DRAG  CURVE 

must  be  input  using  cards  ji» 


CARD  H 

-  OMIT  WHEN  I T APE= 1  • 

COLUMN 

I  TEM 

SCALING 

1-5 

NGPSl J) 

=  NUMBER  OF  CLASSES  OF 

FRAGMENTATION  DATA  IN 

15 

POLAR  ZONE  J.  ALL  ZONES  MUST  HAVE  SAME  NUM3ER 

OF  CLASSES.  (MAX  OF 

10) 

6-10 

BLANK 

15 

11-20 

THSU ( J ) 

=  UPPER  ANGLE  DEFINING 

POLAR  ZONE  J  FOR  WARHEAD 

FIO.O 

(  DEG) 

21-30 

THSL ( J) 

=  LOWER  ANGLE  DEFINING 

POLAR  ZONE  J  FOR  WARHEAD 

FIO.O 

(  DEG) 

CARO  I  -  OMIT  WHEN  ITAPE=1. 


COLUMN  ITEM  SCALING 

1-10  FM(J.K)  =  FRAGMENT  MASS  FOR  POLAR  ZONE  J,  AND  CLASS  K  F10.0 

(GRAIN) 

11-20  F  v ( J » K )  -  INITIAL  VELOCITY  FOR  POLAR  ZONE  J  AND  CLASS  K  F10.0 

! FT/SEC) 

21-30  F F  A ( J  *  K ) =  FRAGMENT  PRESENTED  AREA  FOR  POLAR  ZONE  J  AND  F10.0 

CLASS  K  <SQ.  IN.) 

3 i -AO  blank  fio.o 

41-50  RHO (  J  »  K )  =  FRAGMENT  DENSITY  OR  NUMBER  OF  FRAGMENTS  FIO.O 

51-60  D1(K)  =  WIDTH  OF  RECTANGULAR  FRAGMENTS  IN  CLASS  K  (IN)  FIO.O 

61-70  D2(K)  =  HEIGHT  OF  RECTANGULAR  FRAGMENTS  IN  CLASS  KUN)  FIO.O 

Dina  AND  D2(K)  are  required  only  if  fragment  mass  is 
GEING  OPTIMIZED  or  PARAMETERIZED 


CARD  1  IS  REPEATED  NGPS(I)  TIMES. 

CARDS  H  AND  I  ARE  REPEATED  FOR  EACH  POLAR  ZONE 

CARD  II  (  OMIT  IF  MINI  .NE.  1  > 


COLUMN  ITEM  SCALING 

1-10  FM1  =  FRAGMENT  MASS  FOR  MINI  PLANE  OF  FRAGMENTS.  FIO.O 

(GRAINS) 

11-20  F V 1  °  INITIAL  VELOCITY  FOR  MINI  PLANE  OF  FRAGMENTS.  FIO.O 

I  FT /SEC ) 

21-30  FPA1  *  FRAGMENT  PRESENTED  AREA  FOR  MINI  PLANE  OF  FIO.O 

FRAGMENTS.  (SO.  IN.) 

31-40  BLANK 

41-50  RHO 1  «  NO.  OF  FRAGMENTS  IN  MIN)  PLANE  OF  FRAGMENTS.  FIO.O 

51-60  DD1  “  WIDTH  OF  RECTANGULAR  FRAGMENTS  IN  MINI  PLANE  FIO.O 

OF  FRAGMENTS.,  (  1  N»  ) 

61-70  Dt'**  ■  HEIGHT  OF  RECTANGULAR  FRAGMENTS  IN  MINI  PLANE  FIO.O 

OF  FRAGMENTS.  UN.) 


IF  INOCO  IS  LESS  THAN  3*  OMIT  CARDS  J1  (SPECIAL  DRAG  TABLE) 


C^RD  J1 

COLUMN  ITEM  SCALING 

ENABLES  INPUT  OF  FLECHETTE  DRAG  COEFFICIENTS 
1-10  CD! I)  =  DRAG  COEFFICIENT  AT  VELOCITY  =  0  FT/SEC  F10.0 

11-20  CD  I  I )  =  DRAG  COEFFICIENT  AT  VELOCITY  =  447  FT/SEC  FlO.O 


<  • 

61-70  :D(7)  =  DRAG  COEFFICIENT  AT  VELOCITY  =  1005  FT/SEC  FlO.O 

CARD  J1  IS  REPEATED  UNTIL  25  DRAG  COEFFICIENTS 
HAVE  BEEN  INPUT(7  PER  CARD).  CORRESPONDING  VELOCITIES 
MUST  BE  0.447,558, 670, 782,893. 1005,1116,1228, 1340, 

1451,1563, 1786, 1898. 20 10 ,2233 ,26 79, 30 14, 3350 ,3908, 
4679.5024,5582,6475,10048. 
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SECTION  I  I B 


COLUMN 

1-72 


COLUMN 

5 

6-10 

11-20 

21-30 

35 

41-50 

51-60 

61-70 

71 


72 


(WARHEADS  WITH  VARYING  CHARGE  TO  MASS  RATIO  - 
OPTIMIZATION  OR  PARAMETERS  RUNS  FOR  SPHERICAL 
FRAGMENT  MASS  ONLY) 


CARD  E 


FORMAT 


ITEM 

WTITLE  =  WARHEAD  TITLE 


CARD  F 


FORMAT 


ITEM 

IWHDT  =  WARHEAD  TYPE*  1  =  SPHERE  OR  TABULAR  DATA* 

2  =  CYLINDER 

NPOLAR  =  NUMBER  OF  POLAR  ZONES  REQUIRED  TO  DESCRIBE 
FRAGMENTATION  DATA  FOR  WARHEAD  (MAX  36) 

PBEXP  =  BARE  CHARGE  OF  TNT  EQUIVALENT  EXPLOSIVE 
IN  THE  WARHEAD  (LBS) 

RHOMTL  =  DENSITY  OF  FRAGMENT  MATERIAL.  BLANK  IF 
USING  STEEL  (LBS. /CU.  IN.) 

MINI  =  1*  INDICATES  WH  CONTAINS  A  PLANE  OF  FRAGMENTS. 

WLGTH  =  LENGTH  OF  WARHEAD  (INCHES) 

WDIAM  =  DIAMETER  OF  WARHEAD  (INCHES) 

FACTK  =  K  FACTOR  ( FT*«2/GR**2/3 ) 

I  FRAG  »  FRAGMENT  DENSITY  INDICATOR  =0.  THE  FRAGMENT 
DENSITY.  RHO( I  * J )  IS  INPUT  IN  FRAGMENTS/ 
STERADI AN.  IF  1  FRAG  «  1*  FRAGMENT  DENSITY 
IS  INPUT  AS  THE  TOTAL  NUMBER  OF  FRAGMENTS 
AND  FRAGMENT /STERADI AN  IS  COMPUTED. 

INDCD  3  DRAG  INDICATOR*  1  =  SPHERE*  2  =  CUBE.  3  - 
RANDOM.  IF  INDCD  GT  3*  SPECIAL  DRAG  CURVE 
MUST  BE  INPUT  USING  CARDS  II 

NOTE  -  IF  I  RAO  »  1*  PBEXP  NOT  REQUIRED 

NOTE  -  LEAVE  FACTK  BLANK  IF  FRAGMENT  PRESENTED  AREA 
IS  USED 


SCALING 

12A6 


SCALING 

I  1 

15 

F10.0 

F10.0 

II 

F10.0 

F10.0 

F1Q.0 

II 


1 1 
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CARD  G 


CARD  G  IS  REPEATED  NPOLAR  TIMES 


COLUMN 

Htn 

1-5  NGPSCI)  =  NUMBER  OF  CLASSES  OF  FRAGMENTATION  DATA  IN 
POLAR  ZONE  I.  ALL  ZONES  MUST  HAVE  SAME 
NUMBER  OF  CLASSES*  MAX  OF  10 

6-10  BLANK 

11-20  THSU(I)  =  UPPER  ANGLE  DEFINING  POLAR  ZONE  1  FOR 
WARHEAD.  MAX  OF  180 ( DEGREES) 

21-30  THSU(I)  =  LOWER  ANGLE  DEFINING  POLAR  ZONE  FOR  WARHEAD 
MAX  OF  180 ( DEGREES ) 


SCALING 

15 


FlO.O 

F10.0 


COLUMN 

5 

6-10 

11-20 
21-30 
3  1-40 

41-50 


CARD  HI 

ITEM 

IFORM  =  0*  TABULAR  WARHEAD  DATA*  1  =  GURNEY  FORMULA 
NMAS  =  NUMBER  OF  MASSES  IN  TABLE  MAX  OF  14, 

NOT  REQUIRED  IF  IFORM  =  l 

ENERGY  =  ENERGY  CONSTANT*  NOT  REQUIRED  IF  IFORM  =  n 
PACKEF  =  PACKING  EFFICIENCY  EP  ,NOT  USED  IF  1FORM*0 
RHOMTX  =  DENSITY  OF  MATRIX  MATERIAL  (LBS/IN**3) 

NOT  REQUIRED  IF  IFORM  =  0 
RHOEXP  =  DENSITY  OF  EXPLOSIVE,  NOT  REQUIRED  IF 
IFORM  =  0 

CARDS  H2-H4  NOT  REQUIRED  IF  IFORM  =  1 


SCALING 

II 

15 

FlO.O 

FlO.O 

FlO.O 

FlO.O 


CARD  H2 

COLUMN 

1-10  SFMAS(I) 


ITEM 

FIRST  MASS  IN  TABLE ( 3RAI NS ) 


SCALING 

FlO.O 


CARD  H3 

COLUMN  ITSM 

1-10  _  SVEL(I)  =  SPHERICAL  FRAGMENT  VELOCITY  IN  CLASS  I  SCF 

( FT/SEC !  MAX  OF  SEVEN  CLASSES 

CARD  H4 

1-10  SRnOU)  -  NUMBER  OF  SPHERICAL  FRAGMENTS  IN  FIRST  F 

CLASS  MAX  OF  SEVEN  CLASSES 

REPEAT  CARDS  H3  AND  H4  FOR  EACH  POLAR  ZONE 
REPEAT  CARDS  H2  THROUGH  H4  NMAS  TIMES 

IFINDCD  IS  LESS  THAN  3.  OMIT  CARDS  U  (SPECIAL  DRAG  TABLE) 
CARD  11  (SAME  AS  CAROJ2  IN  SECTION  IIA) 


scaling 

FlO.O 


FlO.O 
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SECTION  IIC 


(KDGN=1»  FOR  DESIGN  MODE) 


CARD  E 

column 

ITEM 

SCALING 

1-72 

WTITLE  =  WARHEAD  TITLE 

1  2A6 

CARD  i 

COLUMN 

ITEM 

SCALING 

1-10 

RHOEX 

DENSITY  OF  EXPLOSIVE  (LBS./CU.  IN.) 

F10.0 

1 1-20 

RHOSTL 

= 

DENSITY  OF  METAL  (LRC- ; 

LF AYE  ZlrXK  IF  MtTAL  IS  STEEL 

•  *  U  •  V 

21-30 

GC 

- 

GURNEY  CONSTANT  (FT. /SEC.) 

F10.0 

31-40 

= 

DETONATION  VELOCITY  (FT. /SEC.) 

FlO.O 

41-50 

THET 

= 

SPRAY  ANGLE  (DEG.) 

FlO.O 

= 

BLANK,  IF  NOT  HELD  CONSTANT 

FlO.O 

55 

NCOMB 

= 

COMBINATION  NUMBER 

15 

NCOMB  TABLE 

( NCOMB  is  the  flag  used  to  determine  which  sft  of 
DESIGN  PARAMETERS  IS  TO  BE  OPTIMIZED  OR  PARAMETERIZED) 


NCOMB  - PARAMETERS 


1 

WWT 

C/M 

L/D 

H/W 

FM(  1 ) 

2 

WWT 

C/M 

DWH 

H/W 

FM  (  1 1 

3 

WWT 

L/D 

DWH 

H/W 

FM  (  1 ) 

4 

WWT 

L/0 

VOL 

H/W 

FM(  1) 

5 

WWT 

DWH 

VOL 

H/W 

FM  (  1 ) 

6 

C/M 

L/D 

DWH 

H/W 

FM(  1! 

7 

C/M 

L/D 

VOL 

H'W 

FM(  .) 

8 

C/M 

DWH 

VOL 

H/W 

"Mil) 

9 

*  VOL 

WW  T 

C/M 

H/W 

FM<  1  > 

10 

*  VOL 

DWH 

L/D 

H/W 

FM(  1) 

NOTE— 1)  IN  USING  THE  DESIGN  MODE  THE  USER  MUST  ACCOMPLISH  THE 
FOLLOWING 

A.  SELECT  A  COMBINATION  NUMBER  (NCOMB  T A  ILE ) 

BASED  ON  WAR*  EAD  CONSTRAINTS, 
o.  OPTIMIZE  OR  ARAMETERIZE  ANY  VARIABLE  IN 
THE  CHOSEN  COMBINATION  BY  USE  OF  CARD  C. 

NOB-  VM.  ALPHA*  AND  ZONEU  CAN  ALSO  BE 
CPTIMIZED  08  PARAMETERIZED. 

C.  any  VARIABLE  In  THE  CHOSEN  COMBINATION  THAT  IS 
N' <T  OPTIMIZEC  OR  PARAMETERIZED  MUST  HAVE  A 
CONSTANT  VALUE  ENTERED  ON  CARD  G. 

2)*COMBINATION  NUMBERS  9  AND  10  HAVE  NO  UNIQUE  SOLUTION 
FOR  THE  TWO  UNKNOWNS.  AND  SO  CANNOT  BE  USED,  THEY  ARE 
LISTED  HERE  FOR  COMPLETION  ONLY. 
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CARO  G 

(INPUT  VALUES  ONLY  FOR  THE  VARIABLES  TO  BE  HELD  CONSTANT) 


COLUMN  ITEM  SCALING 

1-10  WVOL  =  WARHEAD  VOLUME  (CU.  IN.)  F10.0 

11-20  WWT  =  WARHEAD  WEIGHT  (LBS.)  F10.0 

21-30  DWH  =  DIAMETER  OF  WARHEAD (IN.)  F10.0 

31-40  RID  =  RATIO  OF  LENGTH  TO  DIAMETER  OF  WARHEAD  F10.0 

41-50  RCM  =  RATIO  OF  CHARGE  TO  METAL  OF  WARHEAD  F10.0 

51-60  RHW  =  RATIO  OF  HEIGHT  TO  WIDTH  OF  FRAGMENTS  F10.0 

61-70  FM (1*1)  a  FRAGMENT  MASS  (GRAINS)  FlO.O 


NOTE-  IF  A  VARIABLE  MUST  BE  HELD  CONSTANT  IN  A  GIVEN  RUN. 

THAT  VARIABLE  MUST  APPEAR  IN  THE  LIST  CORRESPONDING  TO 
THE  COMBINATION  NUMBER  INPUT. 


SECTION  III 


TARGET  SECTION 


CARD  K2 

COLUMN  ITEM  SCALING 

5  NELMTS  =  NUMBER  OF  TARGET  ELEMENTS  (MAX  OF  5)  I  „ 

6-10  NESH  =  NUMBER  OF  MATERIAL  TARGET  ELEMENTS  WHICH  ARE  SAME  15 

HEIGHT.  THESE  SHOULD  BE  GROUPED  TOGETHER  AS  THE 
FIRST  TARGETS  TO  BE  INPUT.  THOSE  ON  TAPE  SHOULD  BE 
ARRANGED  IN  ASCENDING  ORDER  OF  FILE  NUMBER. 

15-35  LTYPE(I)  FOR  1=1... .*5  =  0  FOR  MATERIEL.  515 

*  1  FOR  PERSONNEL 


CARD  L 

COLUMN 

1  I  TAPE 


2-3  ISKPP 
4-9  TITLE 


CARD  M  -  OMIT  WHEN  ITAPE=1. 

COLUMN  ITEM 

1-72  ETITLE  =  ELEMENT  TITLE 

NOTE.  IF  THE  TARGET  IS  PERSONNEL,  SKIP  TO  CARD  U 

CARD  N  -  OMIT  COL  1-40  WHEN  I  TAPE*  1. 

COLUMN  ITEM  SCALING 

1-5  BLANK 

6-10  NFMS(I)  a  NUMBER  OF  FRAGMENT  MASSES  USED  TO  DESCRIBE  THE  15 
VULNERABILITY  OF  TARGET  ELEMENT  I  (MAX  OF  10) 

11-15  NVELU)  =  NUMBER  OF  FRAGMENT  VELOCITIES  USED  TO  DESCRIBE  15 
THE  VULNERABILITY  OF  TARGET  ELEMENT  I  (MAX  OF  12) 
16-20  NELVSU)  o  NUMBER  0*  ELEVATION  ANGLES  USED  TO  DESCRIBE  THE  15 

VULNERABILITY  OF  TARGET  ELEMENT  I  (MAX  OF  7) 

23  I KT (  l ) =  0.  THE  TARGET  VULNERABILITY  DATA  AS  A  FUNCTION  OF  12 
ELEVATION  ANGLE  IS  USED. 

1,  THE  TARGET  VULNERABILITY  DATA  FOR  THE  UPPER 
HEMISPHERE  IS  USED. 

2.  THIS  OPTION  INCORPORATED  INTO  THE  PROGRAM  TO 
EVALUATE  VULNERABILITY  FOR  FUTURE  ELEMENTS. 


NOT  NEEDED  AT  THIS  TIME. 

31-40  iTin  »  HEIGHT  OF  TARGET  ELEMENT  CENTROID.  F10.0 
41-50  TLGTHUJb  LENGTH  OF  ELEMENTU),  F10.0 
51-60  TWDTH( ! I eWlDTN  OF  ELEMENT! I )•  F10.0 
61-70  THGTUJ  e  HEIGHT  OF  ELEMENT  ( I  >  •  F10.0 


SCALING 

12A6 


ITEM  SCALING 

=  1,  INDICATES  CERTAIN  VARIABLES  WILL  BE  READ  II 

FROM  TAPE. 

=  0,  NO  VARIABLES  READ  FROM  TAPE 

=  FILE  NUMBER  TO  BE  READ  FROM  TAPE  12 

=  ABBREVIATED  NAME  OF  FILE  TO  BE  READ  FROM  TAPE.  A6 


COLUMN 

ITEM 

SCALING 

1-10 

BLSTLl(I)  =  EITHER  THE  IMPULSE  LEVEL  OR  THE  BLAST 
RADIUS  FOR  PKB  =  1.0  OF  ELEMENT  I. 

F10.0 

11-20 

BLSTL2II)  =  EITHER  THE  IMPULSE  LEVEL  OR  THE  BLAST 

RADIUS  AT  WHICH  PKB  BECOMES  EQUAL  TO  0. 

F10.0 

30 

IRAD(I)  =  0  INDICATES  BLAST  LEVEL  IS  USED 
=  1  INDICATES  BLAST  RADII  ARE  USED 

CARD  P  -  OMIT  TYPES  P  THROUGH  T  WHEN  ITAPE  =  1 

11 

COLUMN 

ITEM 

SCALING 

1-70 

AVtIT.JT.KT)  =  VULNERABLE  AREA  FOR  FRAGMENT  MASS  IT.  14F5.0 

FRAGMENT  VELOCITY  JT  AND  FRAGMENT  STRIKING 

ANGLE  KT  (KT  =  l.NELVS) 

CARD  Q 

COLUMN 

I  TEM 

SCALING 

1-40 

FMAS(I.ITT)  =  AN  ASCENDING  ORDERED  TABLE  OF  FRAGMENT 
MASSES  USED  IN  THE  VULNERABILITY  DATA  FOR 
ELEMENT  U  UTT  =  l.NFMS) 

CARD  R 

8FS.0 

COLUMN 

ITEM 

5CALIN6 

1-70 

'EL(l.JTT)  =  VELOCITY  <JTT=  l.NVEU 

CARO  S 

14F5.0 

COLUMN 

I  TEM 

SCALING 

1-30 

ELV(I.KTT)  =  AN  ASCENDING  ORDERED  TABLE  OF  FRAGMENT 
5TRIKING  ANGLES  USED  IN  THE  VULNERABILITY 
PATH  FOR  ELEMENT  I  (KTT  =  l.NELVS) 

CARD  T 

6F5.0 

COLUMN 

I  TEM 

SCALING 

1-30 

VEETDU I.ITT.KTT)  =  MINIMUM  LETHAL  FRAGMENT  VELOCITY 
FOR  TARGET  ELEMENT  1.  FRAGMENT 

ZONE  N.  FRAGMENT  CLASS  M-IM=1.NGRS> 

6F5.0 

CARO  P  IS  REPEATED  «NFM$m  X  NVEUDJ  TIMES. 

CARD  T  IS  REPEATED  NFMSt I >  TIMES 

CARDS  L  THROUGH  T  ARE  REPEATED  FOR  EACH  MATERIEL  TARGET 


IF  TARGET  IS  PERSONNEL ♦  OMIT  CARDS  N  THROUGH  T  AND  USE  THE 
FOLLOWING  CARDS  IN  THEIR  PLACE 

CARD  U 


column  item  scaling 

1-5  I  TROOP (I)  =  1  FOR  PRONE  TROOPS  15 

=  2  FOR  STANDING  TROOPS>  3  =  FOXHOLE 

6-10  8LANK 

11-20  ZT < I J  =  HEIGHT  OF  TARGET  (FT)  F10.0 

21-30  BLSTLKI)  =  EITHER  THE  IMPULSE  LEVEL  OR  THE  BLAST  F10.0 

RADIUS  FOR  PKB  =  1.0  OF  ELEMENT  I. 

31-40  6LSTL2U)  =  EITHER  THE  IMPULSE  LEVEL  OR  THE  BLA^T  F10.0 

RADIUS  AT  WHICH  PKB  BECOMES  EQUAL  TO  C. 

45  I  RAD ( I )  =  0  INDICATES  BLAST  LEVEL  IS  USED  11 

=  1  INDICATES  BLAST  RADII  ARE  USED 


CARD  V  -  OMIT  WHEN  ITAPt  =  1. 


COLUMN  ITEM  SCALING 

1-10  A  =  CONSTANT  DEFINING  CASUALTY  CRITERION  FlO.O 

11-20  B  =  CONSTANT  DEFINING  CASUALTY  CRITERION  FlO.O 

21-30  C  =  CONSTANT  DEFINING  CASUALTY  CRITERION  FlO.O 


REPEAT  CARDS  M.U.  AND  V  FOR  EACH  PERSONNEL  TARGET 


CARD  W1  OMIT  TYPES  W1  THROUGH  2  IF  IOPPAT=0 

COLUMN  ITEM 

1  NLPAR  a  NUMBER  OF  VARIABLES  TO  BE  PARAMETERIZED  IN 
THE  PATTERN  OPTIMIZATION  SECTION. 


SCALING 
1 1 


CARD  W2  REPEAT  NLPAR  TIMES.  OMIT  IF  IOPPAT  -  0 


COLUMN  ITEM  SCALING 

1-6  NAMES ( M )  »  NAME  OF  VARIABLE  TO  BE  PARAMETERIZED.  A6 

11-20  XMNLIMi  «  VARIABLE  MINIMUM  FlO.O 

21-30  XMXL(M)  »  VARIABLE  MAXIMUM  FlO.O 

31-40  XL  t M )  •-  DELTA  X  FlO.O 


PATTERN  OPTIMIZATION  VARIABLE  TABLE 


COLUMN 

1-5 

1-5 

1-5 

1-4 

1-4 

1-5 


I  tem 

N0LTS  ■  NUMBER  OF  BOMBLETS 
TARHt  «  TARGET  HALF  LENGTH* 

TARKW  »  TARGET  HALF  WIDTH. 

SIGW  =  STANDARD  DEVIATION  OF  DEuiVERY  ERROR  IN  THE  RANGE 
DIRECTION. 

S1GL  •=  STANDARD  DEVIATION  OF  DELIVERY  ERROR  IN  TN£ 
DEFLECTION  DIRECTION. 

AREAL  «  LETHAL  AREA.  IF  THIS  VARIABLE  IS  PARAMETER! ZED 
LAPAR  MUST  BE  1. 
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CARD  X 


OMIT  IF  IOPPAT  =  0 


COLUMN 

ITEM 

SCALING 

1-10 

TARHL  = 

TARGET  HALF  LENGTH. 

F10.0 

11-20 

TARHW  = 

TARGET  HALF  WIDTH. 

F10.0 

25 

IDTAR  = 

0.  INDICATES  ELLIPTICAL  TARGET  SHAPE 

11 

1,  INDICATES  RECTANGULAR  TARGET  SHAPE 

CARD  Y 

OMIT  IF  IOPPAT  =  0 

COLUMN 

I  TEM 

SCALING 

1-10 

PATHL  = 

WEAPON  PATTERN  HALF  LENGTH. 

F10.C- 

11-20 

PATHW  = 

WEAPON  PATTERN  HALF  WIDTH. 

F10.0 

25 

I  DP AT  = 

0.  INDICATES  ELLIPTICAL  PATTERN  SHAPE. 

1.  INDICATES  RECTANGULAR  PATTERN  SHAPE. 

I  1 

31-40 

SIGL  * 

STANDARD  DEVIATION  OF  DELIVERY  ERROR  IN  THE 
RANGE  DIRECTION. 

F10.0 

41-50 

SI  GW  = 

STANDARD  DEVIATION  OF  DELIVERY  ERROR  IN  THE 
DEFLECTION  DIRECTION. 

F10.0 

51-60 

DL 

DELTA  ON  PATTERN  SIZE 

F10.0 

61-70 

N8LTS  = 

NO.  OF  BOMBLETS. 

110 

CARD  Z 

COLUMN 

ITEM 

SCALING 

1-10 

FP(l) 

=  WEIGHTING  FACTOR  FOR  TARGET  ELEMENT  I. 

F10.0 

11-20 

• 

FP(  2) 

, 

• 

41-50 

FP  ( 5 ) 

NOTE.  SUM  OF  THE  ABOVE  MUST  BE  1 


NOTE.  RUNS  MAY  BE  STACKED  BY  PLACING  A  *RUN  CARD 

(*  IN  COL.  1*  RUN  IN  COLS.  7-9)  AFTER  LAST  DATA 
CARD  AND  THEN  REPEATING  THE  CARDS  IN  TABLE  I. 

THE  GTITLE  CARD  IS  USED  TO  INDICATE  IF  THE  W H 
DATA  AND/OR  THE  TARGET  DATA  ARE  TO  BE  CHANGED. 

IF  NEW  WARHEAD  DATA  IS  IN  GTITLE  CARD  COLS.  1-16* 

THEN  TABLE  I  I A  OR  I  IB  CARDS  MUST  BE  INPUT. 

IF  NEW  TARGET  DATA  IS  IN  COLS.  20-34,  THEN  THE 
ENTIRE  TABLE  III  CARD  DECK  MUST  BE  INPUT,  COLS.  35-72 
MAv  CONTAIN  ANY  OTHER  INFORMATION  TO  BE  PRINTED  AT 
TOP  OF  PAGE. 

THIS  MAY  BE  DONE  ANY  NUMBER  OF  TIMES,  THUS  ENABLING 
THE  USER  TO  OPTIMIZE  SEVERAL  WARHEADS  AGAINST  VARIOUS 
TARGETS. 


SPECIAL  INPUT  SECTION  USED  TO  STORE  DATA  ON  TAPE 


CARD  1 


%  COLUMN  ITEM  SCALING 

1-6  WT  =  ABBREVIATED  NAME  OF  THE  FILE  TO  BE  READ  FROM  A6 

TAPE. 


CARD  2 

COLUMN  ITEM  SCALING 


1-72 

WTI TLE 

a 

WARHEAD  TITLE 

12A6 

CARD  3 

:olumn 

ITEM 

SCALING 

1-10 

NPOLAR 

5 

SAME  AS  CARD  M  BUT  READ  FROM  TAPE. 

1  f> 

ol-70 

factk 

a 

SAME  AS  CARD  H  BUT  READ  FROM  TAPE. 

F10.0 

71 

1  FRAG 

c 

SAVE  AS  CARD  H  BUT  READ  FROM  TAPE. 

F10.0 

72 

INDCD 

a 

SAME  AS  CARD  H  but  read  from  TAPE. 

FlO.Q 

CARO  4 

COLUMN 

I  rev 

scaling 

1-5 

7  jPSCJI 

* 

NUMBER  OF  CLASSES  OF  FRAGMENTATION  DATA  IN  IS 

POLAR  ZONE  J*  ALL  ZONES  MUST  HAVE  SAME  NUMBER 

OF  classes. 

6-10 

Blank 

11-20 

TMSUt J> 

£ 

UPPER  ANCLE  OfcFtNiNO  POLAR  ZONE  J  FOR 

WARHEAD  r  10.  ft 

IQEOI 

21-30 

ThSL  *  J I 

a 

LOWER  ANCLE  DEFININ'.  POLAR  ZONE  J  f  OR 

WARHEAD  F10.0 

toeoi 


card  $ 


COLUMN 

SYMBOL 

t  tem 

scaling 

1-10 

f,mi  J.K  ) 

m 

same  as 

CARO  M  OUT 

READ  FROM  TAPE. 

FIC.O 

n-20 

PV» J»K> 

a 

SAME  AS 

CARO  M  OUT 

READ  FROM  TAPE. 

F  10*0 

21-30 

AOOM 

a 

blank  - 

F  10.0 

31-40 

XDUM 

t 

Blanks 

F  10.0 

41-50 

Rhoi J.K) 

*  SAME  AS 

i  CARD  H  OUT 

READ  FROM  TAPE. 

FI  0.0 

51-60 

BOUH 

« 

BLANKS 

F  10*0 

61-70 

COUM 

BLANKS 

f  1 0 . 0 

CARD  C 

COLUMN 

1-6 

TAfiGTl 

ITEM 

•  abbreviated  HAKE  OF 

SCALING 

FILE  TO  BE  READ  FROM  TAPE.  M 

si 


CARD  7 


COLUMN 

SYMBOL 

I  TEM 

SCALING 

1-72 

ETITLE  = 

ELEMENT  TITLE 

12A6 

CARD  8 

COLUMN 

SYMBOL 

ITEM 

SCALING 

1-5 

BLANK 

6-10 

NFMSt  I  ) 

=  NUMBER  OF  FRAGMENT  MASSES  USED  TO  DESCRIBE 

THE 

15 

VULNERABILITY  OF  TARGET  ELEMENT  I.  (MAX  OF 

10) 

1 1-15 

NVEL  (  I  ) 

=  NUMBER  OF  FRAGMENT  VELOCITIES  USED  TO  DESCRIBE 

15 

THE  VULNERABILITY  OF  TARGET  ELEMENT  I  (MAX 

OF 

21  ) 

16-20 

NELVS ( I ) 

=  NUMBER  OF  ELEVATION  ANGLES  USED  TO  DESCRIBE  THE  15 

VULNERABILITY  OF  TARGET  ELEMENT  I  (MAX  OF 

10) 

25 

IKK  I  )  = 

0.  THE  TARGET  VULNERABILITY  DATA  AS  A  FUNCTION 

OF  I  1 

ELVS  ELEVATION  ANGLE  ARE  USED. 

= 

1.  THE  TARGET  VULNERABILITY  DATA  FOR  THE  UPPER 

HEMISPHERE  ARE  USED. 

31-40 

ZT<  I  )  = 

HEIGHT  OF  TARGET  ELEMENT  CENTROID  (FT) 

FiO.O 

CARD  9 

COLUMN  SYMBOL  ITEM  SCALING 

1-70  AV(IT.JT.KT)  =  VULNERABLE  AREA  FOR  FRAGMENT  MASS  IT.  14F5.0 

FRAGMENT  VELOCITY  JT  AND  FRAGMENT  STRIKING 
ANGLE  KT  ( K T  =  l.NELVS) 

CARD  10 


COLUMN  SYMBOL  ITEM  SCALING 

1-40  FMAS( I  *  I TT )  «  AN  ASCENDING  ORDERED  TABLE  OF  FRAGMENT  8F5.0 

MASSES  USED  IN  THE  VULNERABILITY  DATA  FOR 
ELEMENT  I.  <  KT  =  l.NFMS) 

CARD  11 

COLUMN  SYMBOL  ITEM  SCALING 

1-70  VEL  (  I  *  JTT )  *  VELOCITY  (J.'T«  A.NVtL)  14F5.0 

CARD  12 

COLUMN  S\  .80L  ITEM  SCALING 

1-30  EL'M  I »KTT )  =  AN  ASCENDING  ORDERED  TABLE  OF  FRAGMENT  oFS.O 

STRIKING  ANGLCS  USED  IN  THE  VULNERABILITY 
PATH  FOR  ELEMENT  I  <  KTT  »  UNELVSI 


CARD  13 

COLUMN  SYMBOL  ITEM  SCALING 

1-30  VEET8L1 1 » ITT  tKTT)  «  MINIMUM  LETHAL  FRAGMENT  VELOCITY  6F5.0 
FOR  TARGET  ELEMENT  I.  ZONE  N.  FRAGMENT  MASS  M. 

ZONE  N>  FRAGMENT  CLASS  M-< M« 1 ,NG*S ) 

CARD  14 


COLUMN  SYMBOL  ITEM  SCALING 

J -10  A  »  CONSTANT  DEFINING  CASUALTY  F1Q.0 

11-20  8  CONSTANT  DEFINING  CASUALTY  FI 0.0 

21-30  C  ■  CONSTANT  DEFINING  CASUALTY  F1Q.0 
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APPENDIX  II 


GURNEY  EQUATION  CORRECTION  FACTOR 

In  developing  the  algorithms  for  this  modification,  it  was  discovered 
that  there  is  a  degree  of  uncertainty  about  the  Gurney  formula.  The  Gurney 
formula,  derived  in  Reference  4,  predicts  the  fragment  fly-off  velocity 
of  a  warhead  as  a  function  of  explosive  characteristics  and  charge-to-metal 
ratio.  The  uncertainty  concerns  the  accuracy  of  the  fly-off  velocity  pre¬ 
diction.  Several  publications  (Reference  2,  for  example)  suggest  using  70 
or  80  percent  of  the  value  predicted  by  the  Gurney  equation.  There  is  some 
disagreement  in  this  uncertainty,  however.  References  5,  6.  and  7  mention 
that  for  a  warhead  with  a  length-to-diameter  ratio  of  less  than  2,  the 
Gurney  formula  yields  high  results  and  as  a  result  should  not  be  adjusted. 

In  an  attempt  to  include  recent  and  accurate  state-of-the-art  technology 
in  this  modification,  an  investigation  was  performed  to  determine  the  most 
recent  and  the  most  accurate  algorithm  to  predict  fragment  fly-off  velocity. 

Reference  5  presents  some  experimental  data  illustrating  the  variation 
of  fragment  fly-off  velocity  as  a  function  of  warhead  length-to-diameter 
ratio.  However,  a  representative  equation  is  not  given;  therefore  the 
results  cannot  be  easily  utilized.  The  purpose  of  this  appendix  is  tc 
derive  such  an  equation. 

The  graphical  data  in  Reference  5  approximate  the  general  form  of  an 
inverse  exponential  curve  approaching  unity.  For  this  reason,  the  following 
general  equation  was  chosen  to  fit  the  data: 

Y  «  l~aebX  (I 1-1) 

where  Y  is  the  dependent  variable  (correction  factor)  and  X  is  the  independent 
variable  (length-to-diameter  ratio).  The  constants  a  and  b  must  be  determined 
and  they  can  be  computed  using  the  least  squares  method. 

Equation  ( 1 1 - 1 )  must  be  transformed  into  a  polynomial  equation  before 
the  least  squares  method  con  be  applied.  This  can  be  done  by  using  the 
following  transformation. 


ln(l-Y)  »  In  a  ♦  b  X  In  s  (II -2) 

ln(i-Y)  «■  lr.  a  ♦  b  X 

Equation  (ST-2)  is  in  the  first  order  polynomial  forte, 

U»  AfBV  (H-3) 

S3 


where  A  and  B  are  constants  and  U  and  V  are  variables.  In  regression 
analysis,  it  is  shorn  that  an  equation  can  be  fit  to  a  set  of  data  points 
in  the  form  of  equation  (II-3)  by  use  of  the  normal  equations, 

£U.  =  nA+BEV. ' 

1  1 


and 


EU.V.  *  AEV.+BEV. 2 

IX  11 

where  n  (i  *  1,  ...  ,  n)  is  the  number  of  data  points,  A  and  B  correspond 
to  the  A  and  B  in  equation  (II-3),  and  (U. ,  V.)  is  the  ith  data  point. 
Solving  for  A  and  B,  11 


n 


n 


(I 1-4) 


nEV.U.-EU.EV. 

B  2  2~ 

i 

By  applying  equations  (II -4)  to  equation  (II-2),  the  constants  in  oquation 
(II -1)  can  be  determined 


EX^nO-Y.HElnO-YpKiX) 


nlX^-CEX^2 


a  «=expl 


/ElnU-Y^-bD^ 


\ 


n 


) 


(11-5) 


A  computer  program  was  written  to  perform  these  calculations.  The 
Gurney  equation  correction  factor,  y*  was  found  to  be 


Y  -  1-.4486  e“l'2345<L/D>  (II-6) 

The  program  source  statements  and  the  output  are  listed  on  the  next  following 
computer  printout  sheets. 
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14/09/70 


GURFAC  FORTRAN  EXTENDED  VERSION  2.0 


PROGRAM  GURFAC (OUTPUT , INPUT,T  APE5=INPUT) 

DIMENSION  X<1Q0),Y<1QO),LINE(1QO) 

1  FORMAT (5F10 .3) 

2  F ORMAT ( 2F10 . 3) 

20  FORMAT (1H1,5X,  *INPUT  DATA (D )  COMPUTED  VALUE(C)*) 

21  FORMAT  (  8X,*X*,9X,  *Y*-,  6X  >*L/Q*  ,6X,*GAMHA  -  GAMMA-Y1*^ - - 

22  FORMAT(iHl,//,6X,*0*,18X,*.2*,18X,*«>4*,i8X,*.6*,iaX,*.8*,17X,*l.*) 

31  F0RMAT(iX,5F10.3)  . . -  - 

35  FORMAT*///, 19X,*C0MPUTED  EQUATION--*, /,19X, ♦GAMMA  =  1.0  -*,F8.5,* 

1EXP(*,F8.5,*  L/0)*» ///////, 19 X,*INPUT  DATA  FROM  REFERENCE  5,  PG  4- 
2183.*) 

50  FORMATtiX»F4-*2»lX,iOOAl>  - - - - - 

65  FORMAT (1H*,5X,100A1) 

READ<5, 1)  YMAX,YMIN,DELX,XMAX,XMIN  -  . . . 

ASSIGN  25  TO  M 

CALL  EOF(M)  . — . 

C  INPUT 

10  REAO (5 , 2)  X (N) , Y (N) 

N  »  N  +1  . .  — . 

GO  TO  10 

C  SUM  PROPER  VARIABLES  - -  -- 

25  SUMLNZ  a  0. 

N  *  N  -1 . -  - . - .  ■  - - - 

SL'MX  =  0. 

SUNX2  =  0.  .  -  - 

SMLNZX  =  0. 

DO  15  I  a  1  ,N 
Z  =  l.-Y(I) 

SUMLNZ  •  SUMLNZ+ALOO«Z)  - - 

SUMX  =  SUMX  +  XU) 

SUMX2  *  SUKX2  *  XU)*X(I) 

15  SMLNZX  =  SMLNZX  ♦  ALOG(Z)*X(I) 

C  COMPUTE  CONSTANTS 

AN  =  FLOAT (N) 

C  »  (AN*SMLNZX-8UMLNZ*5UNX>/(AN»SUNX2-SUNX*SU#X4- 
A  =  EXP  (SUHLNZ/AN  -  C*SUMX/AN) 

C  OUTPUT 

PRINT  20 
PRINT  21 
00  30  I  a  1 »N 

VC  -  1.-  A»eXP<C*XUM  . 

YOIFF  =  YC  -YU) 

30  PRINT  31»X(I),YU),XU)  ,YC»  YOIFF 
PRINT  35, A, C 
C  PRINT  GRAPH 

PRINT  22 

SPAN  IftOWWMAXHMtXMh - 
XXsXHIN 

00  40  I  *  1,100 
40  LINE (I)  **  IN 
42  LINEU)  *  1H* 

00  45  I  »  20,100,20 


14/09/70 


GURFAC  FORTRAN  EXTENDED  VERSION  2,0 


45  LINE (I)  =  1H* 

IX  =  IFIX<SPAN*-<i.-A*£XP{C*XX)-YMIN>  +  .51  - 

IF(IX.GE.l.ANO.IX.LE.lOO)  LINE(IX)  =  1HC 
PRINF  50»XX» (LINE  01) » 1  =  1 » ICQ)  -  ‘ - 

IF  ( IX  . GE. i«  AND. IX  «LE. 100)  LINE<IX>=iH 

oo-w-r-=~tn<r - - -  —  . . — - - 

IF < X ( I)  .GE.XX.ANO.Xa>  .LT.XX+DELX)  GO  TO  60 

55  CONTINUE  - .  . 

GG  TO  75 

68  IX  =  IFIX<SPANMV<n-YMXN»  F.5) 

LINE(IX)  =  1HD 

. PRINT- *'5,  . 

LINE(IX)  =  1H 
GC  TO  55 
75  XX  =  XX+OELX 

IF(XX.tr.XMAX)  GO  TO  42 
CALL  EXIT 

EM) - - - 


INPUT 

X 

OATA (O) 

Y 

COMPUTED 
L/C- — 

VALUE (C) 
- -GAMMA-  — 

GAMMA1- Y" 

.150 

.630 

.150 

-  .627  - 

-.003 

.250 

.580 

.250 

.671 

.091 

.350 

.800 

.350 

.709  - 

-.091 

.500 

.720 

.500 

.758 

.038 

.650 

.854 

.  .650— 

. i'799-  — 

--•*.055- 

.950 

.880 

.950 

.861 

-.019 

1.000 

.707 

1.000 

. .869 

.162 

1.000 

.770 

1.000 

.869 

.099 

1.000 

.824 

1.000- 

.869 

.045 

1.000 

.865 

1.000 

.869 

.004 

1.000 

.892 

■  ■ 

- .069 - 

- .  823 

1.000 

.898 

1.000 

.869 

-.029 

1.250 

.970 

1.250 

■  -  .904 

-.066 

1.500 

.640 

1.500 

.930 

.090 

1.500 

.938 

1.500 

.930 

-.008 

1.500 

.950 

1.500 

,930 

-.020 

COMPUTED  EQUATION-- 

GAMMA  =  1.0  -  .<*4060  EXP<-1. 23449  L/O) 


INPUT  OATA  FROM  REFERENCE  5,  PG  4-183. 
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APPENDIX  III 


DISCUSSION  OF  POLAR  ZONE  ASSUMPTION 

It  is  assumed  in  this  and  most  similar  computer  programs  that  no 
significant  accuracy  is  lost  if  the  upper  polar  zone  angle  (0 ^  in  Figure 
I I I - 1)  is  represented  by  the  angle  computed  in  Shapiro's  formula  (0  in 
Figure  III-l).  If  valid,  this  assumption  will  simplify  the  computations 
performed  within  the  program.  In  this  appendix,  the  conditions  under  which 
this  assumption  is  valid  will  be  discussed. 

Figure  III-l  illustrates  a  cylindrical  warhead  with  the  polar  zone  angle 
'(6.)  and  Shapiro's  angle  (0)  pictured.  For  the  assumption  to  be  valid,  it 
must  be  shown  that  arc  length  B  is  very  much  larger  than  the  sum  of  arc 
lengths  A  and  C.  That  is,  it  must  be  shown  that  assuming  the  fragments 
emanate  from  the  center  of  the  warhead  is  an  accurate  approximation  of  the 
real  case  in  which  the  source  of  the  fragments  is  the  side  of  the  warhead 
and  not  a  point  at  its  center. 


Figure  HI-1.  Polar  Zone  Assumption 
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This  assumption  can  be  stated  as: 
B»B+A+C , 


or 

B»A+C 

Since  ij |l2l |13  and  14| |15 
some  obvious  identities  are 

6=ei=82  A*d 

O<0<j  B=r0  C«~WHL 

Now  equation  (III-l)  can  be  stated  as 


r3»d+|lVHL 

The  length,  d,  can  be  determined  from  the  following 


tan  8,  » 


a 


1  |oWH 
1 


a  =>  2  (DWH)  tan  6 


now 


b  »  -  WHL  -  a 


»  ~  WHL  -  7  (DWH)  tan  0j 


tan  6 


2  c 


so 


.  b 
e  tan02 


i  NHL  -  ~  DWH 
tan  0, 


(HI-1) 


(III-2) 


(III -3) 


(I I I -4) 


(HI-5) 


sin  8  =  — 
c 


d  =  c  sin  8 


j  WHL-(DWH)  tan  8j 
2  tan  82 


sin  8 


(III-6) 


Since  8  =  8^  =  8^ 

,  1  WE  -  (DWH)  tan  8  , 

d  =  —  - -  x.. .  -  sin 

2  tan  82 

=  ~  [WHL  -  (DiW)  tan  9]  cos  8 


(HI-?) 


=  j  [WHL  (cos  8)  -  DWH  (sin  6)] 
Combining  equation  (III-2)  and  (III-7) 


0  r»  j  [WHL  (cos  0)  -  DWH  sin  0  +  WHL] 

»  i  [WHL  (1  +  cos  0)  -  DWH  sin  8] 
Since  0  is  constant  for  a  given  situation. 


(III-8) 


r»  Iq-  [WHL  (l  +  cos  9)  -  DWH  sin  0]  (II 1-9) 

Equation  (III-l)  has  now  been  transformed  into  equation  (I I I -8) -  The  simula¬ 
tion  of  fragments  as  emanating  from  the  center  of  the  warhead  becomes  more 
accurate  as  the  distance  of  the  target  from  the  warhead  and  tho  fragment 
spray  angle  become  larger. 


61 

(The  reverse  of  this  page  is  blank) 


/ 


UNCLASSIFIED 

Security  Classification 


DOCUMENT  CONTROL  DATA  -  R  &  D 

(Security  ctaaaitieatioo  of  title,  body  of  abstract  and  indexing  annotation  must  be  entered  when  the  overatt  report  is  clatsilied) 


\  ORIGINATING  activity  ( Corporate  author) 


T 


2a.  RtPORT  SECURITY  C  L  AS5I F  1C  A  t  ION 


Weapon  Systems  Analysis  Division 
Air  Force  Armament  Laboratory 
Eglin  Air  Force  Base,  Florida  32542 


Unclassified 

*6.  GROUP 

_ KlA _ 


J  REPORT  TITLE 


CYLINDRICAL  WARHEAD  DESIGN  OPTIMIZATION 


4  OESCRIPTIVE  HO TES  (Typo  ot  report  and  inctueire  datea) 

Final  Report  (June  -  August  1971) 

3-  AU  T  MOR(S)  (Firmt  name,  middle  initial,  la  at  name) 

Timothy  C.  Ponder 

Frank  A.  Moescher,  Major,  USAF 


a  REPORT  DATE 

la.  TOTAL  NO.  OF  PAGES 

76,  NO  OF  REFS 

March  1972 

67 

7 

am.  CONTRACT  OR  GRANT  NO. 

%a.  ORIGINATOR**  REPORT  NUMBER!*!  | 

6.  PHOJGC  T  NO- 

2543 

AFATL-TR-72-42 

c. 

96.  OTHER  REPORT  NOlli  (Any  other  number*  that  tmay  be  ** signed  1 

<AJa  report)  1 

<2. 

10  OlltRIHUTION  STATEMENT 

Distribution  limited  to  U.  S.  Government  agencies  only;  this  report  documents 
test  and  evaluation;  distribution  limitation  applied  March  1972.  Other  requests 
for  this  document  must  bo  referred  to  the  Air  Force  Armament  Laboratory  (DLYD) , 
Eglin  Air  Force  Base,  Florida  32542. 


II  SUPPLEMENTARY  MOT  St  |J,  SPONSORING  MILITANT  ACTIVITY 

Air  Force  Armament  Laboratory 

Available  in  DDC  Air  Force  Systems  Command 

Eglin  Air  Force  Base,  Florida  32542  _ 

ts  «r5*vi**<r*  ™" 

This  report  documents  the  Cylindrical  Warhead  Design  Optimization  segment  of  the 
Woepons  Optimization  Techniques  computer  program.  This  segment  enables  the  user 
to  optimize  or  parameterize  the  basic  design  parameters  of  a  theoretical  warhead 
for  a  given  target  or  set  of  targets.  The  warhead  lethality  is  determined  as  a 
function  of  tho  basic  design  parameters:  warhead  weight,  warhead  volume,  warhead 
diameter,  charge-to-metal  ratio,  fragment  mass,  ratio  of  warhead  length  to  diameter 
and  fragment  height-to-width  ratio.  This  segment  can  also  optimize  or  parameterize 
height  of  burst,  terminal  velocity,  impact  angle,  and  fragment  spray  angle. 


DD 


ram 

t  WVW 


1473 


UNCLASSIFIED 


Sec*iftt»  UatcilitttiM 


